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I.  OVERVIEW  OF  CONTRACT  EFFORT  AND  RESULTS 
A.  Program  Overview.  Directions  and  Objectives 

This  program  has  as  its  primary  objective  the  quantitative  exploration 
and  generalization  of  the  core-corona  model  of  imploding  plasma  load  dynamics 
(Refs.  1,  2),  in  close  collaboration  with  radiation  physics  and  modeling  work 
ongoing  at  NRL  (Refs.  3-8). 

As  a  qualitative  suninary,  one  may  say  that  the  core-corona  model  arises 
from  four  physical  considerations,  all  interrelated.  The  first  assumption 
is  that  of  a  sharp  density  fall  off  in  the  outer  regions  of  the  annular 
plasma  load.  The  transition  region  between  the  dense  core  plasma  and  the 
halo  of  corona  plasma  surrounding  it  is  associated  with  a  change  from 
classical  to  anomalous  resistivity,  due  to  the  onset  of  marginally  stable 
microturbulence  in  the  low  density  corona  plasma.  The  modes  comprising  the 
turbulence  are  driven  unstable  by  the  electric  field  (i.e.,  the  generator 
voltage)  which  attempts  to  drive  the  electron  drift  speed  above  a  threshold 
value,  but  the  resulting  current  is  held  to  a  marginally  stable  intensity  by 
the  enhanced  resistivity  which  results  from  the  instability.  The  model 
invokes  current-driven  ion  acoustic  turbulence  for  this  mechanism.  The 
instability  sets  in  once  the  current-induced  drift  speed  exceeds  the  coronal 
sound  speed  (V(j^g>  c^),  and  it  develops  quickly  to  a  nonlinear  staqe  of  ion 
acoustic  turbulence.  The  constraint  of  constant  current  density  across  the 
core-corona  boundary  implies  a  tendency  for  the  outer  plasma  to  heat.  The 
ohmic  dissipation  is  nearly  constant  or  slightly  enhanced  due  to  greater 
resistivity,  but  the  radiative  losses  fall  due  to  the  density  drop.  Such  a 
picture  is  physically  sensible  due  to  the  large  axial  electric  fields  in 
wire  plasma  generators,  which  can  quickly  excite  a  tenuous  outer  plasma  to 
instability  in  the  earliest  phases  of  the  implosion. 
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The  second  assumption  is  that  the  high-density  annular  wire/plasma 
core  stops  a  penetrating  coronal  electron  in  distances  short  compared  to  the 
core  dimensions,  so  that  the  coronal  heating  can  couple  to  the  core  and 
soften  the  implosion.  However,  it  is  important  to  note  that  coronal 
electrons  will  not  tend  to  execute  straight  orbits  into  the  core,  due  to  the 
large  magnetic  fields  in  the  current  carrying  zone  on  the  core  surface. 
Previous  decoupling  energy  estimates  (Ref.  1),  based  on  a  simple  test  particle 
picture,  may  be  too  pessimistic  and  criteria  taking  the  magnetic  fields  into 
account  have  now  been  included. 

A  third  ingredient  in  the  core/corona  equations  is  an  isothermal  corona, 
and  for  this  one  must  assume  the  corona  to  be  sufficiently  limited  in  space 
that  equilibration  to  an  isothermal  pressure  balance  can  be  established 
rapidly  as  the  core  annulus  implodes.  The  corona  continuously  re-establishes 
its  pressure  balance  as  the  transition  interface  moves  inward.  This  requires 
the  coronal  scale  height  ’*q(T2)  to  be  such  that  'imp’ 

Cg  the  sound  speed  at  T2,  the  coronal  electron  temperature,  and  a  time 
scale  characteristic  of  the  implosion  process.  When  disturbances  moving  at 
the  local  sound  speed  can  traverse  the  corona  quickly,  the  pressure  balance 
maintains  a  vanishing  coronal  radial  velocity  field,  and  turbulent  or  classi¬ 
cal  relaxation  processes  can  maintain  an  isothermal  system. 

The  final  constraint  in  the  original  core-corona  model  equations  was 
that  of  quasi -static  heat  balance  between 

ohmic  heating  in  the  corona, 

net  deposition  of  hot  electron  energy  into  the  core,  and 

coronal  radiative  losses 
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Linder  the  present  program,  this  constraint  has  been  relaxed.  The  quasi¬ 
static  requirement  has  now  been  removed  by  generalizing  the  original  equa¬ 
tions,  because  it  was  found  that  during  the  collapse  phase  the  heat  trans¬ 
fer  from  core  to  corona  was  too  rapid  to  allow  the  quasi-static  determina¬ 
tion  of  coronal  temperature.  In  addition,  the  flux  of  energy  from  the  core 
to  the  corona  has  been  included,  because  as  the  collapse  occurs  this  flux 
becomes  important;  and  opacity  functions  have  been  used  to  better  describe 
the  radiative  losses.  Corona  radiative  loss  is  included  although  it  is 
usually  a  small  effect  on  the  energy  balance.  The  detailed  physics  of  these 
upgrades  and  their  improvement  of  the  quantitative  modeling  of  the  system 
are  discussed  in  the  Sections  III,  IV,  and  V  of  this  report. 

The  interrelations  among  these  "^our  assumptions,  as  relaxed  in  the  pre¬ 
sent  work,  represent  time-dependent  tests  of  the  model's  viability.  The 
model  has  shown  some  ability  to  predict  the  collapse  radii  and  temperatures 
with  a  gratifying  order-of-magnitude  accuracy.  With  the  inclusion  of  the 
improved  opacity  model,  very  appropriate  radiation  pulsewidths  and  amplitudes 
are  also  produced,  first  by  the  zero-D  model  and  code  and  now  by  the  1-0  versi 
The  small,  fast  computer  program  also  makes  predictions  concerning  the  time- 
dependent  current  partition  between  core  and  corona. 

Scaling  laws  are  presently  under  development  with  the  assistance  of  this 
code.  These  will  provide  design  optimization  information  and  information  on 
the  limitations  of  such  radiation  sources. 
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B.  Relation  to  Past  Work 


The  JAYCOR  theoretical  program  in  support  of  the  NRL  exploding  wire 
effort  has  provided  very  significant  new  directions  in  which  to  seek 
quantitative  explanation  of  the  wire  array  radiation  pulse  and  scaling,  and 
quantitative  models  of  some  of  these  new  phenomena  are  now  being  developed 
and  tested  under  the  program.  In  its  first  stage,  as  an  NRL  subcontract, 
the  program  provided  preliminary  descriptions  of  two  important  new  physical 
processes  which  had  not  been  considered  in  previous  attempts  to  model  the 
radiation  output:  the  siphoning-off  of  current  to  a  low-density  hot  corona 
during  peak  compression,  and  the  possibility  of  diode-region  behavior  at 
necked-off  constrictions  in  the  assembled  plasma.  Both  of  these  new  physical 
elements  have  the  potential  to  significantly  revise  the  earlier  theoretical 
modeling  of  the  radiation  from  the  assembled  plasma. 

In  its  more  recent  stage,  as  a  JAYCOR  program  funded  by  DNA,  the 
emphasis  is  on  exploring  the  coronal  effects,  developing  and  using  analytic 
modeling  and  a  small  computer  code  (an  order  of  magnitude  smaller  and  faster 
than  the  existing  MHO  code,  which  does  not  treat  the  corona).  In  this  stage 
of  the  program,  there  is  an  unprecedented  and  very  successful  close  inter¬ 
action  between  the  NRL  and  JAYCOR  physicists.  This  quantitative  exploration 
of  the  Tidman-Colombant  core-corona  model  has  already  brought  to  light  the 
sensitivity  of  the  model's  output  to  the  time  development  of  coronal  tempera¬ 
ture  and  to  the  opacity  model  used  to  describe  the  core  radiation.  Thus, 
the  plasma  and  radiation  physics  are  being  integrated  in  a  way  that  seems  to 
hold  considerable  promise  for  finally  understanding  the  radiation  pulse  time- 
scale  and  scaling  with  experimental  observables. 
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This  report  describes  the  logic  and  some  of  the  output  of  this  work. 
Starting  from  the  Tidman-Colombant  equations  of  the  core-corona  model 
(Ref.  1),  these  equations  were  generalized  and  solved  iteratively.  When 
this  procedure  frequently  failed  to  give  solutions,  the  equations  were  fur¬ 
ther  generalized  to  account  dynamically  (rather  than  implicitly)  for  the 
time  dependence  of  the  coronal  temperature  and  ionization  state.  This  is 
the  stage  at  which  results  were  reported  in  the  semi-annual  report  (Ref.  2). 
When  this  zero-dimensional  treatment  was  observed  to  give  anomalously  high 
coronal  temperatures  because  mass  flow  into  the  corona  could  not  be  properly 
included  (Ref.  2),  it  was  generalized  further  to  include  a  one-dimensional 
(radial)  grid  of  density  n(r)  in  the  core,  which  was  advanced  in  a  simple 
Lagrangian  fashion  using  a  local  self-similarity  law  discussed  in  Section 
IV  of  this  report.  This  allows  a  consistent  treatment  of  particle  flow  from 
core  to  corona,  which  keeps  coronal  temperatures  at  reasonable  values  and 
prevents  unphysical  explosive  expansion  of  the  corona  to  large  spatial  extent. 

This  important  physical  improvement,  however,  was  achieved  at  the  expense 
of  certain  new  headaches  of  code  design  that  are  inherent  in  dealing  with 
the  localized  forces  and  compressions  essential  to  a  1-D  model.  For  example, 
the  code  now  successfully  damps  round-off  induced  ringing  of  coronal  current 
and  temperature,  and  treats  in  a  simple  but  useful  way  the  transition  from 
locally  self-similar  implosion  to  assembly  on  the  axis.  The  1-D  code  still 
gives  slightly  shorter  higher-power  radiation  pulses  than  the  experiments, 
but  we  believe  that  this  results  from  axial  dependence  oF  the  fields  in  the 
experimental  geometries,  and  consequent  time-delay  of  the  compression  at 
increasing  distances  from  the  cathode,  features  not  shared  by  the  mathematical 
model . 
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As  improved  computation  of  K  and  L  shell  radiative  yields  and  opacities 
have  become  available  from  the  NRL  group  (Refs.  6-8  and  subsequent  work),  we 
have  incorporated  these  via  analytical  fitting  functions  good  to  better  than 
10%. 


Finally,  using  the  explicit  spatial  structure  in  n(r,t),  one  can  calcu¬ 
late  more  precisely  the  expected  radiative  losses  from  the  plasma.  The 
linkage  of  this  model  to  radiation  transfer  codes  presently  operating  at  NRL 
is  then  a  relatively  simple  matter  and  is  already  underway.  Since  our  experi¬ 
ence  has  demonstrated  the  extranely  tight  coupling  between  the  evolution  of 
the  radiative  energy  losses  and  the  evolution  of  the  plasma,  the  use  of  more 
precise  radiative  transport  and  ionization  dynamics  is  crucial  to  any  serious 
evaluation  of  the  model's  utility. 

C.  Principal  Tasks 

This  program  has  as  its  primary  objective  the  quantitative  exploration 
and  generalization  of  the  core-corona  model  of  imploding  plasma-load  dynamics 
(Ref.  1),  in  close  collaboration  with  radiation  physics  and  modeling  work  on¬ 
going  at  NRL. 

Its  principal  tasks  have  been: 

0  Investigation  of  the  conditions  of  validity  of  the  original  quasi- 
adiabatic  core-corona  equilibrium  relations. 

•  Examination  and  modification  of  the  core-corona  energy  transfer 
equation  and  the  genera tor- plasma  coupling  equations. 
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•  Development  of  a  code  using  those  relations  to  test  quasi-equilibrium. 

•  Generalization  to  an  explicitly  time-dependent  code  to  describe 
the  more  general  problem  where  the  original  equations  were  inappli¬ 
cable. 

f  Incorporation  of  the  improved  modeling  of  radiation  and  radiation 
escape. 

•  Further  improvement  and  generalization  to  correct  an  unphysical 
excessive  coronal  heating.  As  implemented,  this  involved  general¬ 
izing  the  code  to  allow  a  radial  density  profile,  advanced  in  a 
"locally  self-similar"  way  on  a  Lagrangian  grid. 

•  Exercising  the  code(s)  to  observe  their  advantages,  limitations, 
and  physical  implications. 

•  Derivation  of  heuristic  scaling  laws  based  on  the  computations 

(yield  and  timescales  vs.  generator  parameters  and  initial  conditions). 

•  Use  of  computational  results  as  guidelines  for  simple  analytic  approx¬ 
imations  and  modeling. 

•  Simultaneous  derivation  of  simplified  analytic  models,  to  allow 
algebraic  scaling  laws  governing  the  approximate  dependence  of  yield 
and  timescales  on  generator  parameters  and  initial  conditions. 

•  Maintaining  close  working  contact  with  both  radiation  computation 
efforts  and  1-D  hydro  code  efforts  (WHYRAC,  WHYRAD)  at  NRL. 

e  Presenting  results  and  cataloging  the  behavior  and  caveats  of  the 


model (s) . 
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Virtually  all  of  these  tasks  stand  accomplished.  We  note,  however, 
that  we  are  still  not  satisfied  that  the  behavior  of  the  codes  represent 
actual  physical  consequences,  and  scrutiny  and  testing  of  the  code  will 
be  continuing  under  the  1980  funding. 

D.  Suinnary  of  Results 

This  section  provides  a  suirmary  of  our  conclusions  from  coordinated 
analytic  and  computational  development  over  the  past  twelve  months.  Where 
indicated,  some  of  these  conclusions  are  tentative  (See  the  following  Sec¬ 
tion  E  on  caveats).  Many  follow  from  the  fundamental  physics  of  the  prob¬ 
lem  and  are  on  rather  more  secure  ground.  But  we  hesitate  to  apply  our  quan¬ 
titative  scaling  laws  based  on  the  code  behavior  to  present  experiments, 
because  of  the  axial  time-dependence  of  the  compressions  observed  in  those 
experiments. 

We  wish  to  highlight  the  following  conclusions: 

9  The  quasistatic  approximation  (Ref.  1)  for  describing  the  coronal 
temperature  tends  to  fail  upon  plasma  assembly.  The  explicitly 
time-dependent  treatment  (Sec.  Ill)  is  more  useful. 

•  A  single  short  radiation  pulse  is  observed,  over  a  timescale  apparently 
equal  to  the  sum  of  the  assembly  time  and  the  radiative  cool¬ 
ing  time  at  peak  density.  The  assembly  timescale  depends  on  the 
thickness  of  the  annular  plasma  core  and  its  peak  collapse  speed. 

•  A  tenuous,  extended  corona  surrounding  the  bulk  of  the  collapsing 
plasma  may  reduce  peak  collapse  velocities  slightly  by  removing 
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some  of  the  JxB  force  to  larger  radii,  and  the  force  reduc¬ 
tion  may  quickly  increase  on  assembly  if  the  core  heats  the 
corona.  This  reduces  the  magnitude  of  the  "notch"  in  the 
total  current. 

•  On  the  other  hand,  a  more  dense  corona  of  smaller  spatial 
extent  may  rob  an  appreciable  fraction  of  the  JxB  force 
throughout  the  implosion,  i.e. ,  the  coronal  current  may  in 
this  case  be  an  appreciable  fraction  of  the  total  current 
at  all  times.  This  can  occur  if  current  is  excluded  from 
most  of  the  core  plasma  by  its  high  conductivity. 

•  Implosion  trajectories  (r,  n,  T,  etc.  as  functions  of  time) 
are  observed  which  either  bounce  and  collapse  (BC),  pause  and 
collapse  (PC)  or  simply  collapse  in  a  monotonic  way  (SC). 

With  the  removal  of  coronal  overheating  that  resulted  from  the 
zero-D  model,  (and  with  the  modified  radiation  package)  the 
simulated  plasma  no  longer  bounces  to  large  radii. 

•  There  are  preliminary  indications  that  as  the  generator  pulse 
timescale  is  shortened,  the  radiative  energy  yield  first  in¬ 
creases,  and  then  decreases  (See  Figure  6).  The  decrease 

may  be  due  to  pre-compression  of  the  plasma  by  the  fast  rising 
field,  before  assembly. 

•  For  fast  collapse  speeds  and  low  initial  temperatures  just 
prior  to  assembly,  plasma  momentum  can  drive  the  density  to 
values  for  which  the  radiated  power  exceeds  the  input  rate  of 
energy  from  assembly.  The  plasma  then  cools  rapidly  once  assembly 
is  complete,  and  need  not  bounce  appreciably. 
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0  Two  competing  effects  influence  the  extent  of  plasma  bounce 
after  collapse:  high  implosion  velocities  tend  to  give  more 
bounce  for  intermediate  peak  densities,  but  they  also  tend  to 
drive  peak  densities  higher,  and  thus  increase  the  abrupt  radi¬ 
ative  energy  loss  that  quenches  the  bounce. 

0  At  the  highest  peak  densities,  the  balance  of  line,  recombination, 
and  Bremsstrahlung  radiation  shifts  toward  Bremsstrahlung  because 
of  excited-state  quenching. 

0  "Refrigerative"  collapse  to  near-solid  densities  and  low  temp¬ 
eratures  is  probably  a  real  physical  phenomena,  insofar  as  a 
second  pulse  of  low-termperature  radiation  often  cools  the  plasma 
until  its  density  reaches  condensation  values. 

0  Since  diode  impedances  usually  hold  up  after  plasma  assembly, 
the  behavior  of  the  diode  current  after  assembly  should  be  quite 
different  depending  on  whether  the  plasma  coalesces  uniformly 
or  in  beads  along  the  axis.  (Such  problems  will  be  considered 
during  the  1980  effort.) 

0  The  calculation  of  particle  flux  from  core  to  corona  was  improved 
so  as  to  be  consistent  with  the  energy  flux  Q^j^. 

As  a  result  of  this  correction,  the  coronal  temperature  T2  remains 
at  more  reasonable  values  (a  few  KeV),  and  the  coronal  current  frac¬ 
tion  1  -  e  tends  to  be  lower  for  given  initial  profiles.  This  in 
turn  acts  in  the  direction  of  shortening  the  collapse  time  and  in- 
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creasing  the  final  collapse  velocities.  Only  larger  and/or 
denser  coronas  can  slow  the  collapse  appreciably  by  robbing  the 
core  of  JxB  force.  With  the  coronal  electron  temperatures  thus 
lowered,  we  do  not  see  the  coronal  blowoff  ("exploding  corona") 
observed  in  earlier  code  runs. 

21  -3 

Collisional  quenching  of  the  line  radiation  at  high  density  (nj~10  cm 

is  now  included  in  an  approximate  way.  The  quenching  at  these  densities  is 

22 

assumed  to  be  so  rapid  that  above  a  critical  density  of  order  10  no  line 
radiation  is  emitted.  Continuum  radiation  is  assumed  blackbody-1 imi ted 
at  the  highest  densities  and  a  continuum  opacity  estimate  is  inserted  at 
intermediate  densities.  Especially  for  the  higher  final  collapse  velocities, 
the  compressed  density  exhibited  by  the  code  can  eventually  reach  the  critical 
value  in  a  large  fraction  of  the  core,  so  that  in  this  approximation  late¬ 
time  line  radiation  is  emitted  only  from  an  outer  layer  of  the  core,  and  in 
lesser  amounts  from  the  (hotter)  corona.  As  the  outer  regions  of  core 
plasma  collapse  onto  this  blackbody  central  region,  they  accrete  as  they 
radiate  their  energy.  This  may  help  explain  the  discrepancy  between  certain 
simple  scaling  laws  based  on  full -volume  radiating  arrays  (Ref.  5)  and 
observed  experimental  trends  (Refs.  6,  7).  Line  radiation  leakage  from 
the  core  with  a  finite  rate  would  then  combine  with  semi-transparent  radiation 
escape  from  the  diminishing  envelope  mass,  with  continuum  radiation,  and 
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with  frequency-dependent  opacity  effects,  to  give  the  final  radiation 
pulse  profile. 

Further  light  will  be  shed  on  this  question  as  the  self-consistent 
radiation/ionization  model  of  Duston  et  al  (under  development)  is  tested 
and  incorporated  into  the  SPLAT  code,  but  the  most  recent  runs  with  the 
present  radiation  package  indicate  that: 

(1)  Continuum  radiation,  because  it  escapes,  contributes  an 
appreciable  part  of  the  radiation. 

(2)  The  rise  time  of  the  radiation  pulse  is  governed  by  the 
implosion  velocity  upon  collapse. 

We  also  emphasize  that  motional  effects  allowing  greater  escape  of  line 
radiation  have  not  yet  been  included  in  the  SPLAT-ID  Code.  The  inclusion  of 
these  effects,  however,  would  only  increase  the  tendency  of  radiative  collapse 
and  cooling. 

The  tendency  of  the  assembled  plasma  to  reach  high  densities  (ilO'^  cm'"^), 
cool  rapidly  by  radiation,  and  have  negligible  energy  remaining  for  bounce  or 
radiation  thereafter,  provides  a  more  physical  interpretation  of  the  singular 
"refrigerative  collapse"  idealization  (Ref.  18).  The  actual  "collapse" 
may  often  be  delayed  by  colli  si onal  quenching,  which  stores  energy  in 
thermal  energy  and  allows  it  to  be  radiated  away  only  at  the  continuum  rate. 
Motional  broadening  would  tend  to  restore  some  of  the  otherwise-trapped 
line  radiation,  but  probably  cannot  greatly  affect  the  collisional  quench¬ 
ing. 
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E.  Caveats 


We  caution  that  the  modeling  described  in  this  report  is  still  in 
its  developmental  stages  and  is  continuing  under  the  1980  funding,  along¬ 
side  investigations  of  noncylindrical  behavior. 

It  is  also  possible  that  another  physical  mechanism  must  be  found  to 
account  for  the  diversion  of  some  of  the  JxB  work  either  into  plasma 
thermal  energy  rather  than  implosion  kinetic  energy,  or  into  kinetic  energy 
of  the  outer  layers  of  plasma,  delayed  somehow  in  their  arrival  on  axis. 

Since  the  classical  viscous  stress  on  the  fluid  core  can  only  restruc¬ 
ture  the  radial  momentum  density  into  a  preferred  profile  resembling  the 
locally  self-similar  solutions  discussed  in  Section  IV,  the  use  of  classical 
viscous  behavior  to  slow  the  speed  of  implosion  is  a  self-limiting  process. 

In  order  to  directly  store  momentum  imparted  by  surface  forces  one  must 
create  a  steep  gradient  in  the  radial  flow  field  V^(r,t),  This  will  generate 
a  very  unstable  situation  and  may  even  force  multi -streaming  within  the 
fluid,  which  will  in  turn  produce  an  enhancement  of  the  effective  viscosity 
coefficient  and  redistribute  the  momentum  over  the  entire  core  plasma.  Once 
redistributed,  this  momentum  sets  the  implosion  speed  of  the  entire  core 
plasma.  In  calculating  the  implosion  speed  using  the  present  model,  which 
assumes  an  infinitely  fast  redistribution  of  the  momentum  density,  one  has 
already  ancitipated  the  end  result  of  the  dissipative  process  detailed 
above.  The  result  is  that  given  surface  forces  can  never  produce  a  slower 
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core  implosion  (than  that  presently  calculated)  unless  momentum  density  is 
allowed  to  build  up  in  the  outer  layers  of  the  plasma  and  remain  there  on 
a  time  scale  characteristic  of  the  entire  run-in  phase.  Even  a  modest 
classical  viscosity  would  force  a  significant  redistribution  of  this  momen¬ 
tum,  to  say  nothing  of  the  effects  of  multistreaming. 

On  the  other  hand,  if  a  second  degree  of  freedom  were  excited  in  any 
attempt  to  induce  radial  motion,  then  the  momentum  imparted  to  this  additional 

mode  would  represent  momentum  unavailable  to  the  radial  motion, and  the  general 
speed  of  implosion  would  be  reduced  for  any  time  history  of  the  external 
driving  fields.  A  good  possibility  for  such  a  second  degree  of  freedom,  and 
one  which  we  propose  to  explore  through  drift-kinetic  theory,  is  the  excitation 
of  vortex  streets  or  convective  cells  in  the  (r,2)  plane.  Convective 
instabilities  are  known  to  exist  for  field  geometries  similar  to  those  of  Z- 
pinch  plasmas  (Refs,  16,  17),  and  there  is  good  reason  to  believe  that  the 
strongly  magnetized  electron  fluid  in  the  outer  regions  of  t.he  core  could  be 
excited  into  a  vortex  street  layer  by  radial  E-fields  arising  as  electrons 
respond  to  local  JxB  forces  and  drag  the  ions  behind.  Such  a  coupling  of 
radial  accelerations  and  vortex  excitation  would  diminish  the  radial  momen¬ 
tum  transferred  to  the  denser  core  regions  and  slow  the  implosion  directly. 
Subsequent  dissipation  of  the  vortex  motion  would  act  as  a  late  contribu¬ 
tion  to  the  thermal  energy  and  participate  in  lengthening  the  radiaton 
pul se. 
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If  such  motion  is  favored  in  these  systems,  the  impact  on  the  radiation 
pulse  length  could  be  significant.  In  the  present  core/corona  model  the 
pulse  width  is  a  very  sensitive  function  of  the  implosion  velocity.  In 
general,  all  mechanisms  which  tend  to  lower  the  implosion  speed  will  also 
lengthen  the  radiation  pulse,  weaken  the  current  notch,  and  produce  the  peak 
radiation  from  a  cooler,  less  dense  plasma. 

The  1-D  locally  self-similar  code  gives  bounce  dynamics  in  much  better 
qualitative  agreement  with  experiment  than  the  zero-D  code.  At  the  high 
densities  achieved  in  the  simulated  plasma  assembly  after  the  main  radiation 
pulse,  the  continuum  alone  could  (without  opacity)  radiate  away  sufficient 
energy  to  abolish  or  greatly  reduce  the  magnitude  of  the  bounce. 

F.  Relation  to  Other  Analysis 

In  the  past  year,  some  analyses  of  imploding  plasma  radiation  sources 
have  been,  and  are  being  carried  on  independently  at  other  laboratories. 
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Physics  International  (PI),  Maxwell  Laboratories  (MLI),  the  Naval  Research 
Laboratory  (NRL),  Systems  Science,  and  Software  (S  ),  and  Science  Applica¬ 
tions,  Inc.  (SAI)  have  all  produced  relevant  papers  on  various  aspects  of 
the  implosion  and  radiation  modeling.  In  this  section  of  our  overview,  we 
give  brief  comparisons  of  some  of  these  analyses  with  our  own  (which  is  dis¬ 
cussed  more  fully  in  Section  II  through  V  below). 

1 .  Relation  to  PI  Analysis 

The  Theoretical  effort  at  Physics  International  (Ref.  20)  has  involved 
heuristic  scaling  laws  based  on  experimental  results,  optimization  and  scal¬ 
ing  of  the  coupling  of  generator  and  plasma  (without  corona),  and  the  1-D 
Lagrangian  MHO  Code  SM/RAD.  A  brief  and  necessarily  partial  summary  is  as 
follows:  Experimentally,  the  radiative  yields  above  about  1  keV  seem  to 
scale  roughly  with  the  square  of  the  peak  current  (as  expected  since  the 
peak  implosion  speed  scales  roughly  as  JxB),  with  lower  values  for  gas  puff 
loads  than  for  wire  loads  (as  expected  because  peak  densities  should  prob¬ 
ably  be  lower  for  the  gas  puffs). 
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To  calculate  a  simple  scaling  of  radiative  yield,  PI  assumes: 

(i)  a  massive  array  with  inductance  constant  during  the 
current  rise; 

(ii)  that  the  circuit  is  crowbarred  at  peak  current; 

(iii)  that  the  final  kinetic  energy  is  the  difference  between  induct¬ 
ively  stored  energy  at  peak  current,  and  that  at  peak  compression. 

(iv)  that  the  radiated  energy  is  the  sum  of  the  peak  implosion 
kinetic  energy  and  the  inductive  energy  of  the  fields; 
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(v)  that  the  ratio  of  initial  to  final  radii  of  the  array 
is  fixed  independently  of  the  machine  parameters 

(vi)  that  the  voltage  pulse  has  a  simple  parabolic  shape. 

-2  -1 

These  assumptions  give  a  scaling  of  the  radiated  energy  as  V  t  Z.  R(b) 

K  ‘9 

where  b=2uQJln( r^/r^) (Z^t^)"^  and  R(b)  rises  from  about  0.6  to  b  =  0  to 
saturation  at  R  1.3  when  b  >  2;  i.e.,  R(b)  is  roughly  proportional  to 
until  Zg  drops  below  some  critical  value.  Here  t^  is  the  generator  voltage 
rise  time  and  V  is  the  peak  voltage. 

The  1-D  SM/RAD  Code  uses  Braginskii  transport  coefficients  and  flux- 
limited  thermal  conductivity,  but  uses  a  variable  multiplier  on  the  electron- 
ion  collision  rate  (usually  set  at  10  to  50  to  give  reasonable  timescales). 

As  does  our  own  code,  it  starts  from  an  ionized  cylindrically  symmetric  plasma, 
and  does  not  describe  the  early  time  history  of  the  plasma  formation.  The 
arbitrary  multiplier  has  the  effect  of  making  the  resistivity  of  the  entire 
plasma  larger  by  an  order  of  magnitude,  and  allowing  complete  penetration  of 
the  current  into  the  (core)  plasma.  There  is  no  coronal  plasma.  One  difference 
then,  is  that  our  JxB  force  occurs  more  on  the  outside,  while  that  in  the  PI 
1-0  Code  is  more  distributed  through  the  core  plasma. 

The  results  of  the  SM/RAD  Code  runs  for  an  initial  array  radius  of  0.5 
cm  indicate  a  bounce  of  most  of  the  plasma  mass  out  to  about  0.2  cm  after 
assembly,  leaving  low  density  on  axis.  By  contrast,  although  the  JAYCOR 
SPLAT  code  has  a  very  crude  collection  algorithm  for  slowing  the  inter¬ 
penetrating  plasma  streams  on  axis  during  assembly,  the  qualitative  effect 
achieved  is  that  of  assembling  the  plasma  at  the  axis  rather  than  bouncing 
it  off  itself. 
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The  SPLAT  and  SM/RAD  Codes  both  appear  to  give  very  short  assembly 
periods,  on  the  order  of  2ns. 

2.  Relation  to  Maxwell  Analysis 

In  the  work  of  Katzenstein  (Ref.  7),  a  figure  of  merit  is  optimized  over 
the  field  of  parameters.  The  figure  of  merit  is  the  ratio  of  d(KE)/dt  at 
the  end  of  the  run-in  phase  to  the  power  that  the  generator  could  supply  to 
an  impedence-matched  load.  The  numerator  is  just  the  power  dissipated  by 
the  motional  impedence  L  .  The  report  does  not  deal  directly  with  the  magni¬ 
tude  or  time  scale  of  the  radiation  pulse  at  peak  compression,  nor  does  it 
predict  the  dynamics  and  timescale  of  core-plasma  assembly  at  the  axis.  The 
ohmic  resistance  of  the. load  is  neglected  in  comparison  with  the  motional 
resistance  resulting  from  time-varying  inductance  L  during  collapse,  and  the 

internal  enerqy  and  pressure  of  the  plasma  during  run-in  are  neglected  com¬ 
pared  to  its  kinetic  energy.  As  in  the  present  work,  the  plasma  is  assumed 
to  be  cylindrically  symmetric  and  roughly  uniform  axially,  and  the  voltage 
Vg(t)  across  the  plasma  is  taken  as  given. 

Both  the  usual  Z-pinch  geometry  and  a  3-pinch  geometry  were  considered 
in  the  Maxwell  analysis. 
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By  comparison,  the  present  report  focuses  in  more  detail  on  the 
physical  processes  in  the  usual  Z-pinch  mode  and  does  not  consider  the 
0-pinch  mode.  It  concentrates  on  the  effects  of  the  corona  of  plasma 
surrounding  the  dense  plasma  implosion;  robbing  the  core  plasma  of  its 
current  and  its  thermal  energy  near  peak  compression,  allowing  a  longer 
cooler  radiation  pulse  from  the  dense  core  plasma  and  reducing  the  magni¬ 
tude  of  the  dip  in  total  current  at  peak  compression.  The  present  work  is 
not  per  se  an  optimization  of  the  coupling  of  plasma  load  to  generator,  but 
provides  the  computational  basis  for  direct  optimization  of  radiation  out¬ 
put  and  spectrum,  without  the  intermediary  of  a  figure  of  merit 
unrelated  simply  to  radiation  output.  While  the  equations  in  the  JAYCOR 
Code  (developed  under  this  Contract)  are  more  involved  and  numerous  than 
those  of  the  Maxwell  circuit  analysis,  the  code  is  still  much  faster  and 
simpler  than  the  usual  MHO  Codes  required  to  model  the  core  dynamics,  and 
provides  the  coupling  of  the  coronal  plasma  effects,  which  cannot  ily 
be  modeled  by  MHO  Codes. 

3.  Relation  to  Mosher  Energetics  Analysis 

Mosher  (Ref.  3)  has  made  a  useful  analytic  estimate  of  titanium  K-line 
radiation  as  a  function  of  the  kinetic  energy  e|^  of  implosion  (per  unit 
length)  just  prior  to  collapse,  the  average  radius  ^  of  the  imploded  plasma, 
and  the  mass  (per  unit  length)  m  of  the  load.  Radiation  is  assumed  "black- 
body"  at  temperatures  below  a  critical  value  and  optically  thin  above  this 
value.  The  line  radiation  decreases  monotonical ly  with  temperature  above 
this  critical  temperature  (about  50eV).  There  is  no  corona  in  the  model. 
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K-line  radiation,  upon  reaching  2.5keV,  rises  abruptly  as  a  function  of 
implosion  kinetic  energy  above  a  theshold  which  depends  on  the  ratio  of 
ki netic-energy- thermal ization  power  to  peak  radiated  power.  Conclusions 
include:  advantage  of  decreased  load  mass  and  recommendation  of  longer 
pulse  and  implosion  times. 

In  contrast,  our  present  still -preliminary  computational  work  appears 
to  indicate  that  the  time-integrated  radiative  yield  has  a  maximum  as  the 
generator  rise  time  is  varied,  first  rising  with  increasing  rise  time  as 
the  peak  implosion  speed  is  increased,  but  then  falling  once  the  generator 
rise  time  is  fast  enough  to  heat  the  plasma  earlier  in  the  assembly  phase 
and  force  its  expansion  before  as  high  a  density  is  achieved.  The  dynamics 
of  the  plasma  motion  in  our  code  are  purposely  somewhat  oversimplified  and 
probably  do  not  describe  plasma  density  profiles  very  accurately,  especially 
for  very  fast  generator  pulses,  but  those  profiles  are  mostly  unimportant 
to  the  radiatiative  yield.  The  experimentally  observed  radiation  pulses 
are  still  somewhat  longer  than  those  given  by  our  code,  but  that  corresponds 
well  to  the  evidence  of  sequential  compression  of  different  axial  points 
at  different  times. 

4.  Systems,  Science  &  Software 

The  S^  group  has  been  carrying  out  2-D  hydrodynamic  modeling  of  imploding 
wire  arrays  (Ref.  10),  as  a  result  of  their  conclusion  that  1-0  modeling  gives 
unrealistically  high  plasma  densities  and  excessively  short  radiation  pulses. 
The  DELTA  code  is  designed  to  provide  accurate  description  of  an  angular  sec¬ 
tor  of  an  imploding  array,  using  a  triangular  grid  and  moving  with  arbitrary 
velocity.  The  code  thus  combines  features  of  the  Lagrangian  (co-moving),  and 
Eulerian  (fixed  grid)  representations,  and  the  equations  of  motion  are  thus 
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written  with  the  nonlinear  (V  )  terms  expressed  as  convolutions.  An 
artificial  viscosity  is  invoked  to  damp  the  large  oscillations  that  would 
otherwise  occur  unstably  in  the  code. 

Although  the  code  may  be  expensive  to  run  for  times  comparable  with 
the  implosion  time,  it  appears  to  be  highly  useful  in  describing  the 
growth  of  important  azimuthal  asymmetries  in  the  early  phases  of  collapse 
(first  10  nanoseconds).  These  asymmetries  then  can  be  preserved  and 
magnified  during  the  remaining  collapse,  and  are  one  strong  candidate  to 
account  for  the  observed  slower  assembly  times. 

By  contrast,  the  JAYCOR  effort  has  been  devoted  to  finding  whether 
other  mechanisms,  such  as  the  coronal  blowoff  plasma,  could  account  for  the 
observed  slower  assembly  times.  Our  preliminary  indication,  still  clouded 
somewhat  by  code  difficulties,  is  that  other  processes  besides  the  corona 
may  be  necessary,  but  there  is  no  reason  to  assume  anomalous  resistivity  in 
the  bulk  of  the  imploding  plasma. 

3 

Other  2-D  studies  by  S  ,  modeling  the  collapse  of  a  cylindrical  array 
of  thin  filaments,  show  the  development  of  helical  shape  to  the  filaments, 
and  there  is  little  doubt  that  this  mechanism  can  be  active  in  a  large  number 
of  observed  shots.  There  appear  also  to  be  some  shots  not  displaying  obvious 
corkscrew  geometry  at  assembly,  and  it  was  the  JAYCOR  mission  to  explore  a 
possible  model  of  their  behavior,  as  affected  by  low-density  current-carrying 
plasma  blowoff,  but  not  azimuthal  perturbations. 
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5.  SAI/NRL/JAYCOR 


Bloomberg  (SAI),  Colombant  (NRL),  and  Goldstein  (JAYCOR)  (Ref.  11), 
address  (i)  the  early-time  behavior  of  the  plasma  as  the  plasma  current 
channels  are  established  and  the  skin  currents  formed,  (ii)  the  qualita¬ 
tive  prediction  of  double-diode  effects  when  axial  beading  of  the  assembled 
plasma  occurs,  and  (iii)  the  original  core-corona  model. 

We  believe  that  each  of  these  areas  has  a  high  probability  of  being 
important  to  the  modeling  of  the  implosion  and  radiation  dynamics.  The 
early-time  development,  although  not  yet  described  in  adequate  completeness, 
determines  the  'initial'  state  of  core/corona  profile,  current-carrying 
layer  size,  and  thicknesses  of  core  and  coronal  regions.  These  parameters 
then  are  important  in  selecting  among  the  various  classes  of  implosions  and 
radiation  trajectories,  as  discussed  below  in  Section  V.  The  present  JAYCOR 
program,  however,  has  not  dealt  with  the  early  time  history,  but  postulates 
reasonable  outcomes  as  initial  conditions  to  the  modeling  of  the  later-time 
behavior. 

The  formation  of  localized  axial  'diode-like'  gaps,  with  consequent 
axial  acceleration  of  ions  and  electrons,  is  also  a  likely  candidate  to  ex¬ 
plain  several  important  effects,  and  we  will  be  carrying  on  a  more  detailed 
modeling  to  follow-up  on  this  topic  in  the  coming  (1980)  contract  year. 

This  report  covers  the  first  year  of  exploration,  generalization,  and  testing 
of  the  core  corona  model  discussed  in  (iii),  and  originally  formulated  in 
Ref.  1. 
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II.  STARTING  POINT:  QUASI -STATIC  CORE-CORONA  EQUATIONS  (Ref.  1) 


The  full  radial  dynamics  of  imploding  plasma  loads  is  expensive  to 
treat  with  ID  hydro-codes  (such  as  the  NRL  axi -symmetric  code  WHYRAC)  due 
to  the  high  Alfven  and  anomalous  diffusion  speeds  that  occur  in  the  low 
density  corona.  But  a  simpler  model  was  constructed  by  dividing  the  plasma 
into  two  idealized  regions,  namely  the  core  (region  1)  and  corona  (region  2). 

A.  Core  Equations. 

The  high-density  core  is  treated  as  a  spatially  uniform-density  fluid, 
o^(t),  with  classical  plasma  conductivity  and  Joule  dissipation.  The  corona 
surrounds  the  core  and  consists  of  low  density  plasma  driven  by  the  axial 
field  E^  in  a  marginal  micro-instability  state  so  that  the  electron  axial 
drift  velocity  is  approximately  the  sound  speed  in  the  corona.  (The  turbu¬ 
lence  level  in  the  corona  adjusts  to  give  an  electron-wave  scattering  frequency, 
for  which  this  is  true.)  Further,  the  corona  is  assumed  to  adjust  its 
structure  on  a  time  scale  short  compared  to  the  dynamical  evolution  times  of 
the  implosion. 

Although  the  B  field  diffuses  rapidly  through  the  turbulent  corona,  it 
diffuses  relatively  slowly  into  the  classical  core.  The  resulting  partition 
of  total  current  between  core  and  corona  is  determined  self-consistently, 
with  the  result  that  the  coronal  Joule  dissipation  can  be  larger  than  the  core 
Joule  dissipation.  The  coronal  temperature  T2  is, however,  limited  by  micro- 
turbulent  diffusion  (thermal  conduction).  The  radius  a(t)  is  the  boundary 
between  these  regions. 
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In  region  1  (b  ^  r  <_a),  one  has  the  fluid  velocity 


V(r)  =  ( 

r  2pjr  ' 


a^  -  r^] 


where  is  the  core  density,  and  b  is  the  inner  radius  of  the  core  annulus. 
Various  models  for  o^/Sp^  are  possible.  We  determine  this  quantity  by  the 


choice 


This  gives,  in  the  core, 

V(r)  (I.i)  _  ,3^, 

d 

=  A^r"^  +  A^r  (defining  Aj  and  A2)  (3b) 

which  allows  the  annular  shell  thickness  to  expand  at  the  sound  speed  during 
the  run-in  phase,  but  keeps  the  radial  velocity  of  the  inner  edge,  V(b), 
finite  at  collapse  on  axis. 

During  the  run-in  phase  (b  f  0),  the  force  balance  and  energy  equations 
of  the  core  are 


(a^-b")  dt  (  1  b  2  >  2  (32.132)  V  ^  ^  , 


-  -  [i  (a2-r2)  -  r?  in  (^)] 

c2(a2.r2)2  L2  1  1  J 


28 


and 


h  .  3M _ 

TO(a"-r?)  ^re  ^ 


3M  -  M  3  «2  in  a/b  .  .  . 

^r-  =  -  27rp,(aa  -  bo)  —  -^  +  a,  A- 

2m.  at  m.  3t  ^  1  ^2  _(j2  1  2 

*  “  "I  *  »“>]-  (^i  ■  ^i> 

where  a  is  the  core  conductivity 

is  the  core  electron  temperature 
Tj  is  the  core  ion  temperature 
e  is  the  energy  invested  in  ionization  to  degree 
pg,  Pj  are  the  core  electron  and  ion  pressures 
Tg  is  the  electron  cooling  time,  10  A 
M  =  ’ro^(a^-b^)  is  the  total  shell  mass/cm  and 

is  the  total  current  flowing  in  the  core. 

The  constants  and  ^2  are  defined  in  Eq.  (3b).  It  was  assumed  that  the 
total  core  current  is  uniformly  distributed  over  a  skin  depth  a  -  r^. , 
i  .e. , 


Tr{a2  -  r?) 


2Il(r^  -  rj) 
rc(a^  -  rp 


i 


where 


a  :5  r  »  max  of 


(9) 


Once  the  inner  edge  of  the  core  annulus  has  reached  the  axis,  one  has 
b  =  0,  b  =  0,  Aj  *  0,  and  ^2  ~  so  V  »  ar/a.  The  momentum  equation  (4) 
then  becomes 


Ma  =  -  2Tra(P2-Pj)  - 


4Ija 


(a*-r^) 


(10) 


where  now  p(a)  ®  P2  and  p(0)  =  p^. 


3.  Original  Radiation  Emission  Model 

The  radiation  emission  term  ^rad  Eq.  (5)  is  in  general  a  complicated 
function  of  the  plasma  parameters,  geometry,  and  atomic  species.  For  the 
production  of  radiation,  Ref.  1  used  an  expression 


rad 


ia-  -  b-)  (Pg^  + 


P  +  A,p,) 
rec  T  I ' 


where  A.^.  is  a  temperature-dependent  factor  that  accounts  for  reabsorption 
of  the  line  radiation,  P^^.  The  Bremsstrahl ung  and  continuum  recombination 
radiation  Pg^.  and  Ppg^..  together  with  the  average  ion  charge  state  Z,  follow 
from  expressions  used  in  earlier  work.  Only  a  few  dominant  lines  for  the 
plasma  atomic  species  in  question  were  included  in  P^,  e.g.,  averages  over 
K  and  L  shell  lines  for  A]  wires  together  with  the  opacity  factor  A.J.,  were 
provided  by  Davis  and  Jacobs  (Ref.  4). 
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C.  Corona  Equations 

Assuming  Tg  =  T2  »  Tj  the  original  equations  for  the  quasi -static 
corona  were 

(11) 

(12) 

where  the  current  density  j  is  along  the  2  direction.  The  marginal 
instability  current  is 


r 

'2  3r 


-iaB 


4tt  r  3r 


j 


s  e 


and  the  ion  sound  speed  given  by 


(13) 


(U) 


for  ions  of  atomic  weight  A  and  charge  state  Z2(T)  in  the  corona.  These 
pressure  balance  equations  can  then  be  integrated.  Defining  the  symbols 

1 2  =  total  coronal  current  in  amps 
»  total  core  current  in  amps 
T2  *  coronal  electron  temperature  in  KeV 
A,Z2  ■  atomic  weight  and  coronal  ion  charge  state 
a  =  radius  separating  regions  1  and  2 
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Cg  =  3.1  10’(Z2T2/A)^  cm/sec 
Aq  »  2.35  10‘*(T2/nQ)^  . 

the  total  coronal  current,  field  and  electron  density  distributions,  and 
characteristic  scale  height  r^  are  given  after  some  algebra  as 

^  [(21^  -  Ig)2  +  +  (21^  -  Ig)  (15) 


wi  th 


A  *  1.4  10'^  (a  amps 


Ig  =  6.45  10“  T|^gyA/Z2)^  ’ 


(16) 


and 


r8 


^^2 

1 


2  xn 


"  (^) 


n-2 


1  *  (7-) 
0 


(17) 


(13) 


where 


n  =  2  {[3.1  10"^  {tT^  ■  (a^)! 

X  »  0.518  X  10^*  A-nVa^ngZg 

/ T  A 

y^^2  *  6-^5  X  10^ j  J  t-  1  +  n/2) 

\^2  / 

r^  »  a  |x  -  1  -  /x^  -  2x|  »  coronal  scale  height. 
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(19) 


D.  Coronal  Temperature  and  Electron  Deposition  in  the  Core 

The  electron  temperature  required  in  the  above  formulas  must  be  derived 
as  an  approximate  solution  to  the  energy  equation.  Having  chosen  the  corona 
far  enough  out  that  radiation  losses  are  small,  we  have 


F  If  =  jEj  -  -  en^CjEj 


(20) 


where  q^^  is  the  radial  heat  flux.  Integrating  radially  this  gives 


IjEj  *  2™  i  -  Qi2 


(21) 


i.e.,  since  the  corona  is  quasi -static  its  total  Joule  heating  is  balanced 
by  microturbulent  diffusion  of  energetic  electrons  into  the  core.  (It  was 
assumed  that  these  electrons  have  long  mfp's  compared  with  the  core  thickness, 
a-b,  so  that  even  though  B  is  roughly  laminar  in  the  core  they  would  drift 
and  deposit  their  energy  throughout  the  core  region  via  collisions.)  On  this 
basis  the  radial  heat  flux  in  a  corona  of  scale  height  r^  was  then  derived: 


q 


er 


3 

I 


^e^2  '^e'^ew 


(22) 


with  Qg  a  typical  value  of  the  electron  gyrofrequency  eB(r)/nic  in  the  core. 
Thus  an  approximate  expression  for  given  as: 


3  0  Z 


(23) 


(It  turns  out  that  during  run-in  the  corona  is  unmagnetized,  i.e.,  v  >  n., 
but  at  assembly  the  reverse  tends  to  be  true.) 
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It  must  be  noted  that  Eq.  (23)  is  not  really  an  explicit  equation  for 
T2  because  r^  depends  on  n(nQ,a,T2,I^)  and  x(nQ,a,Z2(T2)) ,  and  depends 

on  n^,a,T2  and  Also,  note  in  formulas  (21)  and  (23)  that  is  indepen¬ 
dent  of  r  in  the  corona.  This  follows  from  the  quasi -static  assumption  for 
which  7  X  ^  0. 

The  equation  (Eq.  (22))  has  since  been  replaced  by  a  more  general 
expression,  and  T2,  the  time -dependent  corona  temperature,  is  determined  by 
a  more  general  procedure  than  Eq.  (21). 
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III.  MODIFICATION  AND  EXPLORATION  OF  THE  ZERO-D  CORE-CORONA  MODEL 


In  the  first  quarter  of  the  present  contract  several  modifications  to 
the  original  core-corona  model  were  examined  (Ref.  3)  as  a  means  to  a  more 
physical  treatment  of  the  energy  transfer  between  the  two  regions.  First, 
a  different  result  for  '*'^5  obtained  by  (i)  assuming  the  transition  region 
at  r  =  a  to  be  small  compared  to  an  electron  "mean  free  path"  in  either  the 
classical  or  turbulent  region,  and  (ii)  calculating  the  net  heat  flux  across 
the  boundary  using  a  simple  Krook  model  for  the  relaxation  process  in  either 
region.  The  result  above  (Eq.  23)  is  an  analogy  to  the  classical  heat  flux 
obtained  from  a  temperature  gradient,  but  here  the  corona  is  isothermal  and 
the  correction  to  the  heat  flux  must  come  from  a  treatment  of  the  electron 
density  gradient.  When  the  Krook  model  is  applied  consistently  in  both 
regions,  one  finds. 

q^2  =  (8™^)'^  '^2  I  (24) 

with  K  H  1  -  v^2  (1  ^  ^ew  "2(^)lr  =  a 

3/2 

The  effective  coupling  to  core  parameters  through  the  n^^Tj^  term  was 
previously  neglected  and  represents  the  required  heat  flux  of  cooler  core 

3/ 

electrons  which  drift  into  the  corona.  The  inward  electron  flux,  the  n2(a)T2 
term,  reflects  the  lowest  order  correction  from  the  perturbed  coronal  distri¬ 
bution  function  in  the  presence  of  a  spatial  gradient  in  number  density  and 
a  local  azimuthal  magnetic  field,  Q2- 
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Since  the  coronal  radiation  is  small,  a  useful  first  approximation  was 
to  ignore  radiation  originating  in  the  corona  and  require  quasi-static  energy 
balance,  with  almost  all  the  coronal  ohmic  heating  going  into  the  core  via 
Ql^  5  2'Taq^2- 

E^I  (1  -  3)  =  -  '  2TTa(8TTmg)^  n^Tj^  -  K  02(3)12'  ,  (26) 

A.  Equations  Determining  Coronal  Parameters 

In  terms  of  the  coronal  shape  parameters  n  and  6,  one  can  express  ''2^'''^lr=a 

as  j  1  -  n/2  +  n/(l+6)|.  If  one  chooses  to  measure  Tj,  T2  in  keV,  I  in  MA, 
in  MV/m,  r^.  and  a  in  units  of  lO'  cm,  n^,  n^  in  units  of  10^°  cm  ,  then  ' 

the  Krook  model  correction  factor  K  is 


2.07(Z/A)^  Tp  1  -  n/2  +  n/1+6 

,  I,  S)  =  1  -  - ^ - - 

’  aEj  [1  +  3.84  X  10'*(Z/A)  I^S^T2/a"E| 


The  corona  density  at  the  boundary,  n(a),  is  given  by 


n^  =  Sl/errc^  (a^  -  rp 


□efi ning 


g  =  (r./a)^  [l  -  (r./a)^j 


one  can  consider  Eq.  (26)  an  equation  for  8»  the  current  fraction  in  the 
core.  In  terms  of  g  and  3  we  can  write 


n(g.S)  =  h  +  4  gS  +  2  (1  +  2g2)  6^+4  g6'  +  6"}^  ,  (30) 

1-S"  ’ 
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for  3,  as  5  was  a  constant,  and  then  using  the  3  value  produced  to  re¬ 
evaluate  i:(S)-  Such  a  simple  contraction  technique  will  converge  uniquely 
to  the  lowest  T2  solution  of  (32)  available,  if  a  solution  exists.  A  failure 
of  the  heat  balance  relation  is  signaled  by  the  production  of  negative  or 
complex  3  in  the  process  of  solving  (36). 
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A  quasi-static  code  was  written  to  solve  these  pressure-balanced 
equations,  coupled  with  Eqs.  (4)  -  (10)  for  the  kinematics  and  a  radiation 
package  for  the  core.  But  as  the  code  development  progressed,  it  became 
clear  that  such  an  implicit  scheme  for  the  calculation  of  corona  temperature 
based  on  Eq.  (26)  was  inadequate  because  corona  heat  balance  tended  to  fail 
at  peak  compression.  The  primary  source  of  this  difficulty  is  the  inclusion 
of  the  heat  flux  outward  (from  core  to  corona)  in  addition  to  the  inward 
flux  (from  corona  to  core)  originally  used  by  Tidman  and  Colombant.  In  many 
cases  this  heating  of  the  corona  by  the  core  is  a  significant  contribution 
to  the  overall  energy  transport  as  the  model  plasma  becomes  more  tightly 
pinched.  When  it  is  dominant,  the  corona  has  no  effective  channels  for  heat 
loss  (its  radiation  output  is  small  in  general)  and  a  quasi-static  heat 
balance,  (Z2  12  +  Z2  T2)  »  0,  cannot  be  obtained. 

A  typical  time  evolution  of  the  model  using  the  quasi-static  T2  scheme 
is  shown  in  Figure  1.  The  implosion  is  driven  by  the  slow  ramp  in  generator 
voltage  and  the  time  scale  is  in  units  of  the  assembly  time,  the  time  for  the 
annulus  to  close.  As  the  implosion  proceeds,  Tj,  the  ion  temperature  (taken 
identical  in  core  and  corona),  becomes  somewhat  larger  than  the  core  electron 
temperature,  T^,  and  the  corona  electron  temperature  T2,  peaks  near  14  keV. 
The  current  tends  to  decrease  at  or  after  assembly  due  to  the  strong  impact 
of  the  implosion  velocity  on  the  effective  inductance  of  the  plasma.  The 
most  striking  feature  of  this  scenario  is  the  sharp  rise  in  plasma  density 
that  occurs  after  the  short  prepulse  of  radiation.  It  is  this  steep  density 
rise  that  is  responsible  for  the  loss  of  heat  balance  and,  in  this  case,  also 
indicates  a  core  plasma  headed  for  a  refrigerative  collapse.  The  evolution 
of  the  current  partition  fraction  3  also  indicates  this  collapse  in  that, 
after  a  sharp  dip  at  t/t,  ~1,  the  current  again  accumulates  in  the  core. 


increasing  the  inward  magnetic  force,  while  (and  the  outward  pressure) 
decays . 

Such  behavior  is  not  seen  in  experiments,  and  an  important  test  of  this 
model  is  its  capability  to  avoid  a  refrigerative  collapse.  Anomalous 
resistivity  of  the  corona  can  be  a  physically  plausible  mechanism  for  softening 
the  pinch,  but  only  if  the  time  sequence  of  heating  is  calculated  properly. 

In  a  dynamic  sense,  one  seeks  corona  behavior  which  will  tend  to  rob  the  core 
of  its  current  as  the  implosion  progresses  and  thus  weaken  the  forces  driving 
the  plasma  inward.  However,  once  the  outward  heat  flux  mentioned  above 
enters  the  model,  this  property  tends  to  be  unattainable  in  the  quasi-static 
model  just  when  it  is  needed  most,  because  of  the  constraint  of  quasi -static 
corona  heat  balance.  For  this  reason  the  implicit  — Tg  scheme  has  been  dis¬ 
carded  and  an  explicit  calculation  of  d(Z2T2)/dt  utilized  instead.  Such  a 
modification  has  a  significant  impact  upon  the  mathematical  structure  of  the 
model,  as  discussed  below,  and  upon  its  behavior  as  discussed  in  the  following 
secti on . 

C.  Explicitly  Time-Dependent  Model 

The  time  evolution  of  corona  electron  energy  Z2T2  can  be  calculated  in 
a  manner  analogous  to  that  used  in  treating  this  quantity  within  the  core 
plasma.  If  the  net  heat  flux  from  the  core  into  the  corona  (per  unit  length) 
is  denoted  Qj2  and  the  direct  radiation  originating  from  the  optically  thin 
corona  is  denoted  P,  _  then  one  writes 

U2  -  lE^d  -  8)  +  Q^2  ’  ^Z,r 

as  the  rate  of  change  for  corona  internal  energy  per  unit  length,  when  the  net 
power  loss  or  gain  is  calculated  as  shown  on  the  right  and  distributed  over 

the  entire  corona  plasma.  This  simple  form  presumes  that  E^,  the  average  field 
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within  the  corona  plasma,  is  not  a  function  of  radius,  a  constraint  that  will 
also  be  used  in  the  circuit  relation.  Let  N2  denote  the  total  number  of  ions 
per  unit  length  in  the  corona,  Z2  the  charge  state.  Then  one  has 


It  'V2)  ■  ^2  Ez  *  .  (38) 

[In  the  implicit  -  T2  version  of  the  model,  this  relation  was  replaced  by  the 
constraint  ^  (Z2T2)  =  0.]  The  ingredients  for  the  evaluation  of  Eq.  (38) 
are  determined  by  the  nature  of  the  corona  pressure  balance  —  specified 
uniquely  by  the  parameters  T2,  I,  and  g  =  (r?/a^)  (l-r?/a^)"^. 

Defining  a  =  2  cT2/eCj  (39) 


(with  Cj  the  sound  speed  in  the  corona),  and 
3q  s  co/I  , 

we  consider  S  as  a  funccion  of  g  and  S^.  Then  (S  is  given  by 

(n/Z  +  1  -  22/ 

^  ^  (n/2  +  23/3^  -  1) 


(40) 


(41) 


Here  n  is  the  exponent  determining  the  rate  of  decay  of  density  with 
increasing  r  in  the  corona,  while  5  is  the  parameter  determining  the  offset 
relative  to  r  *  a  in  the  corona  profile.  In  particular,  the  corona  electron 
density  profile  is  given  by 


n(r)  »  n^ 


(1  *  S)‘  (r/a)''-^ 

[3  +  (r/a)'’]^ 


eirc5,2  (a^  -  r?) 


(42) 


Given  Q^2  4^  determined  above,  one  need  only  calculate  N2  and  P^.  2  using 
the  proper  corona  profile  in  order  to  complete  the  evaluation  of  Eq.  (38). 


The  total  number  N2  at  fixed  T2,  I,  and  g  is  simply 


N2  i  27r  Z2 


r 

J  dr  r 


n(r) 


Tra^n. 


1  .e. 


N2  = 


2[l  +6] 
n 


(43) 


The  coronal  radiative  loss  is  more  tedious  but  does  not  make  a  large  contri¬ 
bution.  Assuming  an  optically  thin  corona,  the  volume  emission  rate  over  all 
frequencies  (n,T,Z)  can  be  evaluated  at  (n^^,  T2,  Z2)  and  then  scaled 
downward  by  [n(r)/nj^  as  one  integrates  over  the  isothermal  corona  region. 
This  yields 
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Using  Eq.  (42),  this  can  be  written  as 


^2,r  “  ^rad  ^'^a  ^2  ^2^*  ’ 


where 


"  ^1-2/n 


F(a,  0,  n)  =  27ra2  / 

n  ''1  [5  +  x]' 


dx 


or 


F(a,  5,  n)  =  ^ 


5 
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(44) 


(45) 


(46) 
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using  the  hypergeometric  function  2^1(3  b|Z).  The  present  version  of  the  code 
calculates  the  values  of  N2,  Qj2»  ^2  p  using  Eqs.  (27),  (28), 

(30),  (37),  and  (40  -  46),  and  forms  the  explicit  derivatives  ^  {1^2'^  using 
Eq.  (38).  The  corona  radiation  loss  P,  „  is  small  compared  with  that  of  the 
core.  It  is  only  large  enough  to  make  a  1-5%  change  in  the  energy  bookkeeping 
at  earlier  phases  of  the  implosion.  But  despite  its  forbidding  appearance, 
it  is  reasonably  economical  in  the  code. 


The  first  analytic  impact  of  the  explicit  treatment  in  evolving  T2  is 
a  modification  of  the  previous  circuit  relation  to  include  the  direct  effects 
of  corona  heating  and  magnetic  energy.  In  this  simple  model  (Ref.  3)  the 
circuit  consisted  of  an  effective  generator  impedance,  Zg,  in  series  with  an 
inductive  and  resistive  load,  Lgi+  (I,  Rpj_,  Lpj_)-  Here  Lg  represents  the 
diode  inductance  and  the  plasma  voltage  must  be  calculated  using  the 
core/corona  plasma  model  configuration.  If  one  imposes  a  driving  generator 
voltage  Vg(t),  (such  as  the  fit  to  machine  parameters  due  to  J.  Katzenstein 
(Ref. 13)),  then  the  total  current  I(t)  and  all  other  plasma  variables  must 
satisfy 


Vg(t)  =  I(t)Zg  ^  Lgi  it)  I 


Kr-j) 


(47) 


Here  Faraday's  law  is  used  to  equate  the  line  integral  of  E^  and  the  change 

in  magnetic  flux  through  a  rectangle  r:  (a  s  r  <  r^^),Z:[-  l/Z  <Z<  +  A/2), 

I  being  the  discharge  length.  If  the  inner  radius  is  chosen  at  r  =  a,  then 

the  resistive  contribution  (j  n  t)  can  be  arranged  into  an  effective  resis- 

d  c 


sir] 


tance  R*  = 


c — ^  .  2  3  A  a  ^  transverse  Spitzer  resistiv- 

rp^  C‘a  c 


ity  of  the  classical  core  plasma,  suitably  time  dependent.  This  is  so  because  the 
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current  density  j(r)  is  taken  roughly  independent  of  r  within  the  skin  depth 
of  the  core  plasma,  r^  Within  the  corona,  the  current  density  is  not 

proportional  to  E^.  and  the  overall  plasma  ‘resistance’  Vp^/I{r)  must  include 
a  contribution  from  the  anomalous  resistivity  in  the  corona.  Since  this 
resistivity  must  itself  depend  on  (and  so  on  a  linear  relation 
V  ,  -  IR  is  impossible.  In  short,  the  corona,  while  resistive,  has  been 

pjc 

assumed  to  be  marginally  stable  with  V  .  =  c  Such  a  "drift  speed  limiter" 

u*6  5*4 

behaves  as  a  constant  current  source  rather  than  as  a  resistor  in  the  usual 
sense. 


In  order  to  use  Eq.  (47),  one  must  now  explicitly  calculate  the  induc¬ 
tive  contribution  using  I(r^:  r^^  >  a)  and  collect  all  terms  according  to 
their  coefficients,  i,  I,  etc.  If  the  corona  solution  for  I(rj^)  is  substituted 
directly,  one  obtains 
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dr. 


l(r^) 


=  5n(f) 
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2n  (1  +  6> 


(48) 


approximating  r  »  r  ,  the  largest  effective  coronal  radius.  The  circuit 
w  c 

relation  now  becomes. 


VgCt)  =  IZg  +  LpI  +  R»  I 


^  (I  ^  (-^)  -  I  a  an  (1  + 


6)}. 


(49) 


Expressing  all  derivatives  within  the  brackets  in  terms  of  changes  in  funda¬ 
mental  model  parameters,  a,  6,  g,  and  I,  one  has 


44 


V  (t)  ■  Z  I  +  L_  I  +  I  ^  R 
g  g  D  * 


.1 


+  I 


22 


•-@)l 

*  S  <8-6,/2)) 


T)  -  2  +  (0/0^) 


jk 

c 


2n  (1 


_ Y 1  +  Tf^  ^  C0  -  0«/2)\ 

,  WnWBog;20)  K*  _d6 _ i_J 


11-2+14.  (0/0g) 


(5a) 


+ 

c 


(1  -  2e/e„)  T|'^  §  -  (1  *  (8  -  g^/2)  11~^  ^  )  /  (1  *  8,8/28) 
11-2+14.  (0/0^) 


vlth 


?=2 

d0 


7,-1  ^  »  2 
'•  dq 


3  ^  fs  +  (1  +  2gS  6  +  30^  g  +  g^l 

1-3^  [1  +  l^  g0  +  2  (1  +  2g^  0^+  l4.g0^  +  0**’] 


2 

_ 0  g  A  2  0g  +  0  ) _ 

[1  +  1*  30  +  2  (1  +  2g^)  0^  +  I4.  g0^  +  0^] 


(51) 


and 


(52) 


This  relation  is  implemented  in  the  code,  using  the  i  calculated 
with  the  value  d(Z2T2)/dt  returned  from  the  corona  heat  transport 
relation,  Eq.  (38),  and  other  model  variables  available  with  each  time  step. 
The  value  of  I  is  computed  at  each  time  step  and  used  to  advance  I(t)  as 
Vg(t)  drives  the  entire  system.  The  inclusion  of  a  and  g  terms  in  the 
circuit  relation  had  been  unnecessary  with  the  constraint  (Z2T2)  =  0  imposed. 
The  proper  treatment  of  them  here  points  out  the  complex  character  of  corona 
corrections  to  the  plasma  inductance  (I  coefficient)  and  the  corona  impact 
upon  the  connection  between  energy  exchange  and  core/corona  current  partition 
(a  coefficient). 
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Presently  the  0-D  implosion  code  SPLAT  incorporates  Eqs.  (50)  and  (38) 
to  determine  the  evolution  of  I(t)  and  T2(t).  In  addition,  several  other 
modifications  to  the  original  core  corona  equations  are  included  as  minor 
optimizations,  corrections,  and  explorations  of  possible  anomalously  dissipa¬ 
tive  processes. 

In  the  coronal  pressure  balance,  the  electron  pressure  is  supplemented 
by  the  ion  pressure  n^  Tj — with  Tj  taken  the  same  in  core  and  corona,  and 
nj(r)  =  Z(T2)'^  ng(r)  through  quasi-neutrality.  The  entire  development  of  the 
coronal  quasi-static  state  is  unchanged,  save  only  the  replacement  ot  =  (T2  + 
Z2'Vj)2c/ec52’^'^-  ^^5). 

As  in  Eq.  (36)  in  Section  II,  the  core  motion  is  specified  by  a  particular 
velocity  field 

=  G^(''/r.)  +  G2(^/r)  (53a) 


which  also  arises  naturally  as  that  solution  of  7 •_)/  ~  ^  *  b^)  with 

no  viscous  stress.  Here  r^  is  used  as  a  convenient  scale  radius  and 
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(b  I  a) 
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ft)  (7!?) 


(53b) 

(53c) 


as  in  Eq.  (3)  in  Section  II. 


This  velocity  field  is  assumed  to  retain  its  characteristic  spatial 
structure  with  G2  given  by  time  dependent  boundary  contidions,  a,  b.  The 

constraint  b  =  a  -  c^  allows  one  to  develop  an  equation  of  motion  for  a(t) 
in  particular: 
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with  p  the  core  mass  density,  n^,  the  spatially  homogeneous  core  electron 
density  and  temperature,  T2  the  coronal  electron  temperature,  c^  the  core 
sound  speed,  Cj^  its  time  derivative  due  to  changes  in  (Z^  T^),  and 


C^TI^ 
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/  C^^+6  ] 


r 


the  inward  coronal  pressure  in  terms  of  the  parameters  describing  the 
(implicitly  evolved)  force  balance.  The  spatially  integrated  heat  balance 
equation  completes  the  set  of  coupled  differential  equations  required  to 
evolve  the  core  fluid  explicitly.  For  the  electron  fluid,  one  has 
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the  net  core/corona  electron  heat  exchange,  and 


(57b) 


rad 


=  -(3-  -  b-)  CPg^  +  P 
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^ine^ 
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the  radiative  loss  term,  based  on  approximations  due  to  Colombant  and  Tonon 
(Ref.  15),  Davis  and  Jacobs  (Ref.  4  ),  and  a  simple  line  radiation  opacity 
model  due  to  Apruzese  (Ref.  7  ).  The  viscous  heating  term  is  given  by 


^is 
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with 


n*  =  0.73  Tli  Ti  T, 


ee 


(57e) 


resulting  from  a  straightforward  evaluation  of  the  classical  expression,  cf. 
Braginskii  (Ref.  12).  A  similar  expression  obtains  for  the  ion  fluid, 
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with  evaluated  for  appropriate  ion  coefficients.  In  these  expressions 
Tge  is  the  usual  Braginskii  collisional  relaxation  time,  M  is  the  total  mass 
per  unit  length  in  the  core  plasma;  and  m^,  the  electron  and  proton  masses; 
A  and  2  atomic  number  and  (temperature  dependent)  charge  state;  the 
classical  transverse  resistivity. 

The  viscous  heating  term  for  the  ion  fluid  may  be  enhanced  during  periods 
of  compression  at  hypersonic  speeds  (a/c,  <-1.5)  to  reflect  the  ion-electron 
energy  relaxation  time  Tj  ^  rather  than  the  classical  like-species  relaxation 
time.  The  physical  picture  here  derives  from  the  fact  that  ion-ion  relaxation 
in  these  situations  exhibits  Coulomb  logarithms  of  order  one,  indicating  a 
rather  fundamental  inadequacy  of  classical  transport  theory.  As  such  ’ons 
tend  to  behave  more  like  a  strongly  coupled  fluid  than  a  classical  plasma, 
momentum  transferred  from  more  tenuous,  faster  flowing  regions  into  more 
dense,  slowly  flowing  regions  will  tend  to  be  distributed  among  several  ions 
at  a  time  in  many-body  collisions.  These  ions  can  then  be  considered  as  a 
quasi-neutral  clump,  transporting  the  momentum  over  a  path  length  dictated 
by  the  time  required  for  electrons  to  dissipate  their  ordered  motion.  Hence 
an  ion-electron  collisional  relaxation  time  is  more  appropriate  for  the 


calculation  of  the  viscosity. 


While  a  detailed  calculation  of  such  effects  is  clearly  preferable,  a 
simple  ad-hoc  application  of  the  basic  idea  here  can  be  quite  useful  in 
directing  future  invest’-^ations.  At  present,  when  desired,  the  viscosity 
coefficient  is  smoothly  allowed  to  rise  to  n*  =  i 3/c, 
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becomes  less  than  -  1.5,  indicating  an  "average"  clump  of  about  10  ions.  The 
effect  is  to  heat  the  ion  fluid  heavily,  with  this  thermal  energy  more  slowly 
coupled  to  the  electron  fluid  and  future  radiative  losses.  However,  if  the 
plasma  thermal  relaxation  time  becomes  so  short  that  separate  ion  and 
electron  temperatures  equilibrate  rapidly,  this  unbalanced  heating  of  the  ion 
component  tends  to  cause  rapid  oscillatory  heat  exchange  between  the  electrons 
and  ions.  When  such  effects  attempt  to  dominate  the  numerical  integration, 
the  simpliest  and  most  appropriate  course  of  action  is  to  leave  this  heating 
source  out  pending  a  detailed  calculation. 

The  calculation  of  ^  ^  transformed  to 

a  direct  expression  for  T^  or  T2  using  any  explicit  functional  model  for 
Z(Tg)  appropriate  to  the  plasma  species  being  investigated.  As  the  power 
warms  the  plasma  during  compression,  the  resistive  magnetic  diffusion  of  r^. 
is  smoothly  replaced  by  a  simple  transport  of  r^  using  the  velocity  field  V(R) 
above.  The  simple  scaling  of  the  magnetic  diffusion  process  (Ref.  1  )  used 
originally  is  now  active  only  in  the  early  stages  of  the  implosion. 

D.  Self-consistency  of  the  Model  Assumptions 

Of  the  three  surviving  assumptions  among  those  considered  at  the  outset, 
the  one  most  deserving  of  some  scrutiny  is  that  of  the  quasi -static  corona 
From  the  explicit  solution  for  n2(r),  one  can  derive  a  simple  criterion  for 
the  spatial  extent  of  the  coronal  region, 

r^~  a(10  +  96)^/'^ 

being  the  radius  for  which  the  coronal  density  has  decayed  to  about 

( 10  +  96) n  .  Typically  one  finds  n  —  5  to  7 ,  6  ~  .4  to  .6  at  peak 
100 

compression  (maximum  coronal  temperature  and  spatial  extent)  so  that  this 
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density  is  about  0.06  n,  or  less.  If  one  computes  how  many  times  a  sound 

d 

wave  can  transverse  the  distance  (r^  -  a)  in  the  time  required  for  the 
implosion  dynamics  to  change  a  by  10%,  the  result  is  usually  of  order  10  -  100 
and  of  order  unity  only  for  the  largest  a  values  on  the  implosion  trajectory. 

There  is  good  reason  to  believe  that  these  largest  a  values  are  somewhat 
unphysical  anyway,  so  the  quasi-static  corona  assumption  appears  quite  consistent 
with  the  other  elements  of  the  model  dynamics.  The  physical  picture  implied 
is  that  the  corona  can  both  expand  to  track  the  imploding  core  surface  and 
also  rearrange  its  internal  state  on  a  shorter  time  scale  through  the  propagation 
of  sonic  turbulence,  in  order  to  continuously  "hover"  about  the  marginal  stabi  lity 
condition  and  isothermal  pressure  balance. 

The  assumption  that  higher  energy  coronal  electrons  can  be  stopped  in  a 
short  distance  of  core  plasma  also  remains  valid  throughout  the  implosion 
process  insofar  as  higher  core  densities  occur  simultaneously  with  the  high 
coronal  temperatures  and  produce  a  very  short  classical  mean  free  path  for  the 
electrons.  This  constraint  is  even  more  easily  satisfied  in  the  later  versions 
of  the  model  because  the  coronal  temperatures  are  somewhat  lower  than  those 
previously  calculated.  This  is  due  primarily  to  the  improvements  in  the 
coronal  kinetic  energy  transport  calculation  discussed  above. 

Finally,  one  must  examine  the  consequences  of  invoking  a  precipitous 
density  variation  on  the  core/corona  interface  and  here  the  results  are  not 
as  good.  The  predominant  effect  of  this  discontinuity  in  n^  is  to  enhance  the 
outward  kinetic  energy  flux  of  core  electrons  very  rapidly  near  the  peak 
implosion  and  it  can  produce  (through  heating)  a  corona  so  great  in  spatial 
extent  that  it  requires  all  the  available  current  to  maintain  its  pressure 
balance.  In  any  physical  system  such  a  density  gradient  must  remain  finite 
and  a  smooth  connection  of  core  and  corona  density  will  result.  This  is  one 
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of  two  primary  motivations  for  expanding  the  model  to  take  account  of  the 
spatial  density  variation  in  the  core,  the  second  being  the  extreme  sensiti¬ 
vity  of  the  model  dynamics  to  the  accuracy  of  the  radiative  loss  calculation. 

A  properly  physical  radiative  transport  and  loss  calculation  cannot  be 
obtained  if  the  density  variations  are  excluded.  The  required  modifications 
to  the  model  are  discussed  in  Section  IV  below. 

E.  Explicitly  Time- Dependent  Code  Results  and  Evaluation 

The  numerical  parameters  and  behavior  typical  of  this  zero-O  model  as 
implemented  are  described  in  a  short  confidential  supplement  to  this  report. 

The  case  treated  is  meant  to  correspond,  however  roughly,  to  wire  plasma 
imploded  with  diode  and  driver  parameters  corresponding  to  the  Maxwell 
Laboratory  devices.  The  driving  voltage  is  that  Vg(t)  due  to  Katzenstein 
(Ref. 13)  offset  by  a  short  delay  time  to  allow  for  the  annular  plasma  formation 
and  multiplied  by  arbitrary  input  factors  (nominally  unity)  to  allow  varying 
the  peak  voltage  and  rise  time.  This  calculation  used  Vg(t)  *  '^Katz^^  ^  ^^nsec' 
The  evolution  of  {a(t),  b(t),  T^{t),  Vg(t),  I(t),  n^(t),  I.|(t)},  i.e., 

the  inner  and  outer  radii,  core  and  corona  electron  temperature,  total  current 
and  driving  voltage,  core  electron  density  and  core  current,  respectively,  is 
shown  in  Figure  (2a,  b).  Here  a{t),  b(t)  are  in  cm,  and  T^  is  in  kev 
n.j,  I,  Vg  are  in  arbitrary  units. 

The  example  shown  here  is  nearly  one  of  a  class  of  system  trajectories  one 
might  call  the  "exploded  corona".  As  the  implosion  proceeds  the  core  current 
fraction  S  begins  to  drop  rapidly  upon  assembly  of  the  annulus.  With  the 
current  then  shifted  into  the  corona  and  the  reversal  of  produce  corona 

heating  rather  than  core  heating,  the  corona  temperature  rapidly  increased.  This 
in  turn  can  further  shift  8  toward  zero  and  the  0-D  model  often  becomes  unstable 
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If  there  is  sufficient  energy  to  be  removed  from  the  core,  this  tendency  to 
drop  3  win  proceed  until  the  corona  is  forced  to  be  essentially  infinite 
in  extent.  An  extremely  expanded  corona  is  not  particularly  plausible 
because  of  the  constraint  that  a  sonic  disturbance  should  traverse  the  corona 
as  implied  by  the  model  trajectory  in  times  short  compared  to  the  required 
hydrodynamic  time  step.  In  an  ideal  case,  one  would  find  the  fraction  s  to 
be  well  above  0.05  over  the  entire  implosion  process. 

The  second  general  class  of  systems  trajectories  is  quite  the  opposite 
case,  which  one  might  call  the  "refrigerator".  In  this  second  scenario,  the 
fraction  6  is  driven  toward  one  as  assembly  occurs.  This  enhances  the  forces 
compressing  the  core  plasma  and  simultaneously  drives  the  coronal  temperature 
downward.  The  plasma  is  driven  to  even  smaller  radii  and  higher  densities 
where  it  radiates  copiously  and  begins  to  cool.  In  this  final  phase,  T^^,  Tj 
and  T2  are  driven  toward  zero,  dropping  by  an  order  of  magnitude  or  more  in 
a  few  nanoseconds. 

In  both  cases  described  above  the  current  trace  is  quite  similar,  the 

decay  of  I(t)  after  assembly  in  Figure  2b  is  due  to  the  time-dependent  nature 

of  the  plasma  inductance.  If  the  implosion  speed  is  very  large  one  finds  that 
2Jta 

the  term,  -  ,  in  Eq.  (50)  can  easily  be  dominant  over  R*  and  Z  .  Any 

c^a  3 

experimental  current  trace  not  exhibiting  this  downward  trend  would  be  inter¬ 
preted,  in  terms  of  this  model,  as  one  corresponding  to  a  slow  implosion 
velocity.  In  any  case,  the  behavior  of  the  circuit  variables  is  not  particularly 
enlightening  as  to  the  reason  for  the  "refrigerator"  and  "exploded  corona" 
bifurcation  among  the  system  trajectories. 
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The  strength  of  the  core  radiation  loss  and  the  coronal  heating  are  the 
primary  variables  controlling  the  transition  from  "refrigerator"  to  "exploded 
corona".  For  any  model  conditions,  a  change  of  the  order  of  1"  in  radiation 
volume  emission  rate  P^g^(n,T,Z)  can  force  the  system  from  one  case  to 
another.  For  larger  emission  rates  the  "refrigerator"  model  occurs,  and,  as 
the  emission  rate  is  smoothly  lowered  (for  all  time  steps),  there  will  exist 
a  critical  value  for  this  quantity  such  that  the  "exploded  corona"  mode 
is  established  for  all  lower  values  of  ,T,Z) .  Such  manipulation  of 

Prad  *1^5  no  direct  physical  significance,  of  course.  On  the  other  hand,  it 
is  of  primary  consequence  in  the  model  evolution  at  the  onset  of  the  radiation 
pulse.  The  nature  and  time  history  of  this  pulse  and  its  dependence  on  some 
of  the  system  or  code  parameters  is  discussed  briefly  in  the  confidential 
supplement. 

An  additional  factor  which  figures  heavily  in  producing  a  near  "exploded 
corona"  trajectory  is  the  initial  radius  of  the  annulus.  In  an  obvious  way, 
arrays  of  larger  radii  will  require  the  thermal! ration  of  more  kinetic  energy 
at  assembly  and  thus  heat  both  core  (and  corona)  to  a  greater  degree.  This 
in  turn  tends  to  drive  5  downward  and  there  will  exist,  for  a  specified  driving 
voltage  waveform  and  radiative  loss,  a  sufficiently  large  initial  radius  such 
that  the  corona  will  always  "explode"  near  the  time  of  initial  assembly. 

These  effects  are  indicative  of  several  things. 

First,  the  larger  implosion  speeds  and  flow  kinetic  energies  are  derived 
primarily  from  an  intrinsic  bias  in  any  spatially  uniform  fluid  description 
toward  the  most  rapid  possible  response  to  the  impressed  JxB.  forces  in  the 
surface  layer.  The  fact  that  these  larger  speeds  are  limited  as  well  as  they 
are  by  the  coronal  current  channel  is  a  measure  of  its  rather  wide  dynamic 
range  in  pinch-softening  capabilities;  but  these  capabilities  are  clearly 
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bounded  by  the  upper  limit  to  8^  forced  by  the  logical  structure  of  this 
sort  of  corona. 

Second,  the  radiation  pulse  length  is  tightly  linked  to  the  overall 
time  scale  of  compression  and  expansion  in  the  core  fluid.  While  better 
radiation  loss  calculations  can  significantly  change  this  pulse  length, 
the  core  density  and  temperature  trajectory  is  the  primary  function  to 

concentrate  upon.  It  is  the  actual  input  to  any  radiative  loss  calcu¬ 
lation  and  is,  therefore,  the  ultimate  determinant  of  the  radiation  pulse 
length.  Any  process  that  bears  on  the  density  and  temperature  trajectory 
calculation  has  a  significant  impact  upon  the  loss  calculation  and  attendant 
radiation  pulse  length,  but  not  conversely.  For  example,  the  details  of  the 
total  current  evolution,  I(t),  will  not  be  particularly  sensitive  to  the  loss 
calculation  (even  though  3(t)  is  quite  sensitive).  On  the  other  hand,  this 
total  current  is  quite  central  to  the  density  and  temperature  trajectory 
calculation.  The  code  replicates  this  conclusion  (cf.  Figure  2b)  with  the 
weakened  notch  of  I(t)  (relative  to  that  in  I.|(t))  through  the  complex  radia¬ 
tion  pulse  generated  by  nj(t),  T.|(t)  during  the  imolosion. 

This  is  not  to  say  that  almost  any  radiative  loss  calculation  will 
suffice,  for  clearly  it  will  not.  One  must  take  into  account  the  effects 
dependent  on  density  and  temp^'rature  profiles  and  take  care  to  evolve  the 
plasma  and  the  radiation  self-consistently.  The  density  and  temperature 
trajectory  calculation  never  decouples  from  the  radiative  loss  calculation 
on  any  time  domain  of  interest  through  the  implosion. 

Because  of  this  tight  linkage  between  the  radiative  pulse  length  and 
density/temperature  trajectories,  the  modifications  discussed  in  the  next 
Section  IV  were  undertaken  as  a  prelude  to  more  advanced  calculations  of 
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radiative  losses.  The  improvements  derived  from  this  more  physical  and 
detailed  picture  help  to  bring  the  model  in  closer  agreement  with  experiment 
by  lowering  the  overall  implosion  speeds  and  by  making  possible  more  accurate 
radiative  loss  calculations,  although  even  further  changes  will  be  seen  to  be 
required  in  order  to  confidently  duplicate  experimentally  observed  results. 

A  third  primary  facet  of  the  model  dynamics  is  the  oscillatory  behavior 
after  the  initial  compression.  This  is  aggravated  by  the  fact  that  there  is 
no  dissipative  mechanism,  within  the  oresent  model,  which  can  sap  the  inward 
momentum  of  the  collapsing  plasma  directly. 

The  plasma  now  behaves  as  a  radiatively  damped,  and  quite  non-linear 
oscillator  ^ — using  the  pressures  built  up  through  thermal i zed  flow  kinetic 
energy,  ohmic  heating,  and  core/corcna  heat  exchange  to  reverse  its  radial 
motion. 

In  fact,  the  core  heating  is  the  only  retarding  influence  on  the  compres¬ 
sion  of  the  pinch,  and  this  is  why  larger  radiative  losses  can  bring  on  a 
refrigerati ve  collapse.  In  a  simple  mechanical  analogy,  the  radiative  loss 
rate  eats  away  at  the  maximum  available  potential  barrier  (the  pressure  built 
up  at  the  end  of  implosion)  so  that  a  given  plasma  with  specified  kinetic 
energy  input  can  find  either  a  positive  or  negative  total  energy  upon  reaching 
the  axis,  depending  on  the  radiative  loss  rate.  The  oscillation  is  heavily 
damped  and  is  not  readily  apparent  experimentally,  but  neither  is  the 
refrigerative  collapse.  The  simplest  alternative  is  some  motion  within  the 
corona  (vorticity  is  a  good  candidate)  that  can  act  as  a  momentum  reservoir, 
fed  during  compression  and  dissipated  afterward.  This  possibility  is  being 
investigated,  but  it  can  be  examined  dynamically  only  once  the  improvements 
of  Section  IV  are  implemented  successfully. 
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IV.  SELF-SIMILAR  ID  CORE-CORONA  CALCULATION 


A.  Mathematical  Model  for  Self-Similar  Core  Compression 
1 .  Basic  Requirements 

The  velocity  field  currently  used,  Eq.  53(a),  is  a  special  case  of  the 
following  simple  theorem  for  compressible  fluid  mechanics. 

Thm  1 .  Given  a  set  of  nested  surfaces  v(_x.t)  =  parametrized  by  spa¬ 
tially  constant  numbers  y.,  let  the  quantity  s  s  -  ^  ■  denote  the 

1  n  *  j 

logarithmic  convective  derivative  of  the  volume  enclosed  by  any  two  such 
surfaces,  and  the  values  -u^-n.y.n  the  normal  velocities  of  such  surfaces  as 

D 

they  move  with  a  fluid.  If  the  flow  field  V.(x,t)  appearing  in  ^  is  now 
decomposed  as 

v(x,t)  =  I(x,t)  +  n  (x,t) 

with  2  X  T  50,  7  Eo;  and  7  is  given  as  any  solution  of  7  x  T  =  7u)  (admitting 

only  rigid  rotor  vorticity),  then  (i)  the  field  y_(x,t)  is  determined  uniquely 

by  the  solution  of 

7>  T  =  -  s  ,  interior  to  v.  . 

^  *  J  » 

and  (ii)  all  such  fields  V  exhibit  no  viscous  stress,  viz.  7*  y  e  o,  whenever 

s  =  s  (t).  The  proof  follows  from  a  simple  application  of  Helmholtz's  theorem 
n  n 

and  the  recognition  of  an  equivalent  Poisson  problem  for  T  with  Neuman  boundary 
conditions  on  a  closed,  simply  connected  surface.  A  second  result  is  immediate 
from  the  continuity  relation; 

Thm  2.  Given  the  flow  field  V()^,t)  above,  the  number  density  exhibits  a 
self- similar  evolution  of  the  following  form: 
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i 


n(x,t)  =  n(x  (x,t),  t  ), 


whenever  is  spatially  homogeneous,  Here  one  has 

S(t.t^)  =  -J  dt^  ^  an  =&n 


and 


w 

^(5.i)  =  Z  +  J  Y^CxCtj^).  t 


;  which  are 


(59) 


(59a) 


(59b) 


equivalent  to  the  differential  equation 
^  n(j<,t)  =  s^(t)  n(^,t). 

Of  course,  it  is  the  equation  of  momentum  transfer  which  determines  the 
time  evolution  of  the  flow  field  and  maintaining  the  specialized  form 

n 

of  V(x,t)  demands  a  specialized  form  for  y_.  Seeking  the  most  general 
transformation  preserving  the  property  that  ±'  'L  ~  ~  i'inds, 


or 


with 


and 


^  V(x,t)  =  v(x,t)  +  a^(t) 


(vCViS.O.t,)  + 


=  exp  J  dtj  Oj  (t^)  , 


(60a  ) 


(60b  ) 


(60c  ) 


(60d  ) 
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These  relations  are  equivalent  to  the  condition  •  X)  =  cii(t),  which  says 
that  the  time  development  of  the  local  compression  v-V  is  independent  of  posi¬ 
tion. 


These  results  can  be  related  for  an  isothermal  plasma  by  requiring  the 
preferred  form  of  ^  ^  to  be  produced  from  the  local  combination  of  pressure 
gradient  and  J.  x  ^  stresses  actually  existing  in  the  plasma.  Neglecting  the 
viscous  stress  because  ^(V-V)  =  0,  consider  the  profile  solving 


J  X  B 


men 


(61) 


~  -1  ~  -1 

with  m  =  m^  +  Z  mj ,  T  =  Tg  +  Z  Tj ,  and  ^  given  as  any  solution  of  7"^  = 

■  ^n^^o^  (Z  ^  X  =  0)  geometry  of  interest.  This  profile  now  generates 

the  preferred  form  of  required  in  Eq.  (60a),  and,  so  long  as  this  struc¬ 
ture  for  ^  V  is  maintained,  the  profile  can  be  evolved  using  Eq.  (59)  -  7.V 
remains  spatially  homogeneous. 


If  external  forces  are  involved,  it  is  not  clear  that  any  such  preserva¬ 
tion  of  V  is  possible.  Let  A,  i  4  and  restrict  attention  to  those  forces 
satisfying  a  relation  of  the  form  discussed  below, 

it  is  possible  to  maintain  a  homogeneous  V-V^  e  •  the  expense  of  some  restric¬ 
tion  on  the  time  dependence  of  Expressing  ^  in  terms  of  its  initial  value, 

with  S^(t,tg)  =  y*  dt-j  S^(ti), 
t_ 


and  the  constraint  on  ^  becomes 
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^0  ^  Dt  -  ■  ■  Ot  ^  ^ 

as  ^  and  A,  evolve  in  time.  Can  any  such  A,  exist? 
n  -iJ  — j 

A  special  case  for  the  profile  ^(x,t)  of  interest  here  is  any  form  j^(_x,t) 
which  is  itself  locked  to  the  fluid  flow,  i.e.,  i(^,t)  such  that 

0(x,t)  =  Sj(t)  J(x,t)  (62) 

from  which  it  follows  through  the  "generalized  Ohm's  Law"  (Ref.  14)  that 


[1  +  S.]  J  =  i  +  (VxB)  c'V 


(  53; 


Further,  for  an  isothermal,  highly  collisional  plasma  (S.(t)v- 
a  nearly  homogeneous  resistivity  n*  =  ^Sj(t)),  one  has 

7x  ~  n*(7xJ)  i  2xE_  +  c'"*  IxVxB 

7xJ  :  -  3,  and  7  x  E  =  -  c  "*3^  B,  thus 

-g*  g  ^-1(3^3  +7xVxB);  (64) 

and  using  the  continuity  relation  for  n  through  ^•)/  “  ■  qY  ^'.2i>‘') »  this  yields 


^  a'B^n  g^(i) 

4- 


(65) 


as  the  implied  constraint  on  ^/n  for  self-similar  evolution.  In  the 

limit  of  large  conductivity  this  relation  forces  B  and  n  to  scale  identically 
with  the  local  compression  of  electron  density,  thus  ^  results  in 


^  (JxB/n)  Sj(JxB/")- 


(66) 
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With  the  results  of  Eq.  (66),  and  the  fact  that 


one  can  derive  (for  a  slab  geometry)  a  simple  constraint  on  the  time  develop¬ 
ment  of  If  a  plasma  slab  is  allowed  to  expand  away  from  its  density 
maximum  on  a  plane  x*,  and  exhibits  an  initial  density  profile  solving  Eq.  (61) 
consistent  with  a  linear  spatial  dependence  of  the  velocity  field  V  =  G.j(t) 

(x  -  x*),  then  self-similarity  can  be  maintained  exactly  whenever 


S;\(t,to) 


(TZ)  V 


(Vo^ 


exp  -  J^dt-j  G^(t.|)  . 


Here  the  point  x*  is  an  invariant  or  fixed  point  of  the  mapping  x°  -*  x(t) 
generated  by  )^(x,t).  This  relation  implies  a  particular  decay  rate  for  the 
current  density  in  the  plasma.  The  scaling  parameter  =  Sj(ti)  must  be 
such  that  it  compensates  for  the  temperature  decay  while  the  ratio  £/n  remains 
constant  on  any  expanding  density  contour  (_x  -  )<*)(t).  The  partial  time 
derivative  of  Eq.  (68)  provides  the  requirement 


vJ  7T  1 


and  this  admits  a  simple  physical  interpretation.  The  isothermal  nature  of 
this  expansion  is  preserved  only  because  the  pdV  work  done  to  move  the  plasma 
is  spatially  homogeneous,  i.e.,  ^'*'e  ~  -  V-V,  ^  ^Tj  s  -  V-V.  This 

implies  the  constraint  Sj  =  -  2(7* V)  +  G-j  =  -(V-V)  and  hence  makes  the  imbedded 
^  compatible  with  the  self-similar  expansion.  Now,  if  B/n  is  constant  on 
~  'i*)(T)  and  Sj  =  -(7-V),  the  time  decay  of  Aj  is  due  solely  to  the  density 
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decay  in  ^  =  n  e  ^  throughout  the  current  carrying  layer,  for  any  current 
profile  J(x.,t).  The  self-similarity  is  preserved  because  Z  ^ 

must  change  in  step  with  n & as  the  drift  speed  ^  is  held  constant  at  its 
original  local  value,  in  order  to  preserve  self-similarity.  Clearly  an  arbi¬ 
trary  time  dependence  for  is  not  compatible  with  the  self-similar  expan¬ 
sion. 


The  self-similar  motions  generated  by  initial  profiles  solving  Eq.  (61) 
can  be  extended  to  curvilinear  coordinate  systems  with  more  modifications. 

In  cylindrical  coordinates  the  only  true  self-similar  expansions  or  contractions 
must  be  centered  on  the  origin;  however,  a  thin  annulus  can  maintain  such  self¬ 
similarity  approximately  while  its  ratio  of  thickness  to  mean  radius  is  small. 

In  general  the  profiles  generateo  by  Eq.  (61)  are  appropriate  to  physical 
systems  which  begin  in  a  very  viscous  state  before  transition  to  plasma  and 
any  period  of  strong  viscosity  in  the  later  evolution  will  always  tend  to 
drive  the  velocity  field  back  to  a  profile  such  that  7-V  is  spatially  homo¬ 
geneous.  This,  in  turn,  tends  to  force  the  density  profile  back  to  its  self¬ 
similar  evolution. 

When  the  core  plasma  is  coupled  to  the  marginally-stable  pressure-balance 

corona,  the  regions  of  self-similar  evolution  are  those  for  which  J  is  approximately 

J..(x  t  )-l/°u/k'^L.  and  the  viscous  heating  is  either  spatially  homogeneous  or 
-0  -0  0  -ab  ab 

negligible.  Moreover,  the  constraint  ~  -(7-y^)  cannot  be  maintained  easily 
because  the  factor  must  respond  to  the  external  equations  that  provide  the 
boundary  value  of  £  on  the  surface  of  the  core. 

But,  since  self-similar  solutions  provide  quite  economical  computational 
models,  it  is  quite  useful  to  let  Eq.  (61)  determine  the  initial  profile 
n(_x,t)  over  the  plasma,  and  approximate  its  evolution  to  remain  self-similar 
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(in  the  sense  of  Eq.  (59)),  with  continuity  required  of  n,  V,  and  T  at  the 
core/corona  interface  defi ned  by  the  constraint  that 


J(x,t) 

'  en(x,t) 


(69) 


This  provides  the  opportunity  for  a  more  physical  calculation  of  the  core/ 
corona  electron  heat  exchange  that  is  so  important  in  the  energetics. 


2.  Application  to  Strictly  Radial  Geometry 

Specializing  Eq.  (61)  to  radial  variations  only  is  relatively  straight¬ 
forward  once  a  model  profile  for  J^(r,t)  is  chosen.  The  uniform  current  density 
being  the  preferred  asymptotic  final  state  of  a  magnetic  diffusion  process, 
we  can  set 


E’ 


E  +  c"^  V  xB  = 
-*  -r  -e 


(1  -  Yj,>0 

0  <  0 
R 


(70) 


with  Yp  =  r/rj  -  1,  and  r^  defining  the  innermost  penetration  depth  of  the 
fields.  Then,  using  the  velocity  field  of  the  original  model  (a  rotation- 
free  special  case  of  Thm  1)  with  a,  b  denoting  the  normal  velocities  of  the 
nested  cylindrical  surfaces, the  solution  of  Eq.  61  can  be  expressed  as 


ni(xR] 


=  {■ 
0  t  ' 


-  H(xr) 


-  U 


-1 


(71) 


in  the  limit  that  p-®  in  Eq.  70.  Here  for  convenience  we  have  let 


F_(x,  u,  o)  :  (1  x)’  exp  ['j(cx  +  x‘/2)] 


(72) 
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and 


U,(C.  u.  ct)  =  F^(x,  u,  o)  dx 

0 

This  solution  arises  from  casting  Eq.  61  in  dimensions! ess  form;  n^  is  a 
reference  mean  density,  and  applying  the  integrating  factor 

exp  -  I  rj  j  dR  V^(R)  ]  , 

where 

m  =  m^  +  Z  ^  , 


and 

R  =  r/r,  using  r,  as  a  convenient  scale  radius. 

i  1  1 

In  Eq.  71  ,  the  function  H(x|^)  is  the  unit  step  function,  of  argument 


^R  ^ 


_  r/r.-l 
1 


The  various  parameters  of  the  model  self-similar  profile  remaining  in  Eq.  71 
are  defined  as 


^a  ^  ^a(«a  ' 

Uo  =  m  rja^  G^/T 
i  1  + 

~  m  f*  j  ct  2  ^2  /T 
;i  H  (e;  rj/n*c)Vn^  T 


(73) 
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r 


and  ,  G2  are  given  by  Eq.  53(b,c,).  The  constant  of  integration  is  deter 
mined  by  the  radial  equivalent  of  Eq.  69,  i.e.,  one  defines  the  core/corona 
interface  by  the  electron  drift  speed  in  response  to  the  local  E'  field.  One 
finds  the  appropriate  to  be 


C 


0 


0’ 


(75) 


with  H  E^/en*(T.|  )c.|  (T^ ).  A  further  requirement  is  that  the  density 
profile  should  normalize  to  envelop  a  specified  total  mass.  In  the  case  of 
Eq.  71  with  fixed  by  Eq.  75,  the  normalization  condition  arising  from  a 
total  number,  Ny,  of  atoms  in  the  initial  plasma  becomes 

Z(T^)Ny  =  2t  r:  (R^-1)  f  dx^  n-jix,^)  [l  +  ,  (76) 


with  Xjj  =  (b/rj-1 )/(a/rj-l )<  0.  This  can  be  reduced  to  an  equation  for  the 
parameter  (R,-l)  that  is  easily  solved  numerically,  i.e. 

a 


[l-(Ra-l)x] 


(77) 


The  denominator  on  the  right  is  an  increasing  function  of  (R.-l)  since  as 

a 

(a/rj-l)  is  increased_,the  absolute  value  of  the  profile  n.|CX|^)  is  forced  to 
intersect  n-j  =  n^  at  a  larger  value  of  argument  of  the  decaying  exponential 
appearing  in  Eq.  (71).  One  need  only  increase  R^-l  for  fixed  (b,rj)  (and 
other  model  parameters)  until  an  intersection  of  f(a)  =  R,-l  and  g(a)  =  RHS  of 

a 

Eq.  (77)  is  achieved. 
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The  density  profile  can  be  evolved  in  an  economical  way  only  by 
following  a,  b,  a,  b,  a,  b  in  a  Lagrangian  sense  in  order  to  keep  track  of  the 
time  evolution  of  all  the  profile  parameters  ,  a^,  ,  and  G2.  This  is 

because  the  fundamental  similarity  variable  is  defined  through  an  integral 
equation, 

.t. 

t"o(r.t)  =  r  +  r°  dt‘V^(r(t‘),  f) 
t 

so  that  one  must  always  follow  the  bounding  zones  of  the  self-similar  solution 
region,  and  any  interior  zones  of  interest,  e.g,  rj(t),  in  a  Lagrangian  sense. 

The  primary  computational  economy  arises  from  the  fact  that  all  time  deriva¬ 
tives  required  to  follow  the  mesh  are  simple  expressions,  derived  from  the  pro¬ 
file  parameters  that  can  be  computed  precisely  rather  than  estimated  from  dif¬ 
ference  techniques. 

The  self-similar  profile  derived  above  will  tend  naturally  to  expand 
away  from  its  density  maximum  and  to  fill  the  available  space  unless  external 
forces  change  its  motion.  In  the  case  of  an  annular  profile  responding  to  a 
specific  model  generator  voltage  trace,  there  are  several  factors  which  alter 
the  self-similarity.  These  can  be  accommodated,  however,  by  simple  models 
for  the  run-in,  assembly  on  axis,  and  later  phases  of  the  implosion.  The 
primary  modification  is  that  due  to  changes  in  the  magnitude  and  distribution 
of  external  force  densities  within  the  plasma.  The  original  imbedded  profile 
of  J  x  B  force  must  respond  to  complex  changes  in  the  E,,  ^  surface  fields 
and  also  reflect  an  ambipolar  radial  electric  field  and  the  associated  pE^ 
stress.  The  final  stresses  moving  the  plasma  annulus  inward  will  be  a  complicated 
superposition:  -Vp  +  pE^  iz  ^  overall  drift  inward  is  developed. 
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the  compression  will  tend  to  be  greater  on  the  interior  surface  (V-V  receives 

-2 

a  weak  contribution  ~  r  )  and  as  the  plasma  annulus  closes  down  on  the  axis 
the  density  gradient  will  reverse  (cf.  Eq.  (67))  and  slow  the  inner  surface 
more  violently.  This  violent  deceleration  will  then  propagate  outward  and 
force  a  conversion  of  flow  energy  to  thermal  energy.  This  "collection"  pro¬ 
cess  onto  the  axis  is  the  second  source  of  alterations  to  self-similar  motions 
in  an  expanding  annulus. 

In  contrast  to  these  mechanisms  altering  self-similarity,  there  is 
the  action  of  viscosity,  which  tends  to  restore  the  homogeneity  of  V- V_ when 
included  in  the  momentum  transport  equation.  In  most  applications  of  classi¬ 
cal  fluid  plasma  theory,  viscosity  is  a  relatively  weak  process  and  is  often 
neglected.  The  basis  for  this  approximation  is  a  collision  calculation  that 
assumes  very  weakly  correlated  plasma  particles,  as  measured  by  the  so-called 
Coulomb  logarithm,  A.  As  a  simple  calculation  will  attest,  A  for  ion/ion 
collisions  in  an  A1  plasma  with  Tj  —  8  eV  -  200  eV,  nj  —  10*®  -  10^®  is  not 
large  enough  to  insure  the  correctness  of  classical  transport  theory.  Except 
for  the  higher  temperature  points  along  the  typical  implosion  trajectory, 
this  parameter  is  less  than  or  on  the  order  of  1.  One  should  therefore  treat 
a  strongly  coupled  plasma,  with  transport  coefficients  that  include  the 
influence  of  many-body  collisions.  The  proper  transport  coefficients  are  not 
available,  and  it  is  necessary  to  model  the  general  effects  of  such  strong 
coupling.  The  model  adopted  thus  far  is  based  on  the  heuristic  assumption  that 
the  stronger  collisional  coupling  does  not  change  the  character  of  the  dissi¬ 
pative  fluid  processes  (heat  conduction,  viscosity,  friction)  but  simply  makes 
them  more  dominant.  The  immediate  implication  is  that  time-asymptotic  config¬ 
urations  for  dissipative  fluid  evolution  are  achieved  quasi-statically  in  the 
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strongly  coupled  core  fluid.  For  internal  energy  changes,  this  implies 
maintenance  of  a  nearly  isothermal  temperature  profile;  for  velocity  field 
changes,  a  nearly  homogeneous  7-^  profile  as  required  for  self-similarity;  and, 
for  the  relative  flow  of  ions  and  electrons,  an  equal  radial  drift  velocity 
for  each  species. 

One  may  treat  the  run-in  phase  by  assuming  that  all  stress  developed 
at  the  outer  surface  a(t),  c"^{J,  x  B„),,  is  transmitted  inward  to  the  location 

“•Z  ~~v  a 

of  the  density  maximum  by  a  self-consistent  combination  (pE^  +  (J^  x  Bg)c”^)  of 
stresses.  This  drives  the  density  maximum  to  accelerate  slowly  inward  with 
some  speed  f^(t)  and  (interior  and  exterior  to  the  plasma  density  maximum  in 
the  frame  moving  at  f^)  the  pressure  gradients  and  viscous  contributions  act 
to  spread  the  profile  with  a  V^{r,t)  exhibiting  a  homogeneous  divergence.  This 
simple  model  for  the  evolving  velocity  field  can  be  summarized  as 


r_  =  0,B,/m  c  n„ 
m  a  a  m 


b  = 


j  7n 
m  n 


together  with  the  constraint 


(78a) 

(78b) 


(78c) 


arising  from  the  7.V  homogeneity.  Here  n^^,  r^^,  f^,  r^^^  denote  variables  at  the 
density  maximum;  n^^,  b,  b,  b,  those  at  the  inner  surface.  Such  a  model  has  the 
advantages  of  computational  simplicity  and  compatibility  with  the  requirement 

^  V  =  a,  V  +  ^.  This  enables  one  to  advance  a  radial  mesh  of  suitable 
Dt  —  1  —  -0 
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spatial  resolution  for  the  accurate  computation  of  radiative  losses. 


The  second  phase  is  that  of  collecting  the  plasma  onto  the  axis,  and 
here  one  must  be  guided  qualitatively  by  a  very  complex  mix  of  phenomena. 

The  instabilities  available  to  the  system  because  of  strong  density  and  flow 
gradients  will  tend  to  drive  a  significant  level  of  microturbulence  in  the 
central  regions,  and  the  colliding  ion  fluid  components  will  be  strongly 
coupled  coll isionally  as  well.  A  detailed  treatment  of  this  sort  of  collapse 
is  not  within  the  scope  of  classical  fluid  plasma  theory;  but  the  net  effect 
will  be  to  mix  the  fluid,  to  dissolve  the  density  and  flow  gradients,  and 
to  thermal ize  the  energy  previously  tied  up  in  the  radial  flow.  At  present 
a  very  simple  model  of  this  "collection"  process  has  been  implemented  and  it 
roughly  describes  the  time  development  outlined  above.  As  the  inner  surface 
{b(t))  meets  the  axis,  the  number  of  particles  flowing  past  a  small  collection 
radius,  r^;  is  computed;  and  the  flow  energy  of  the  local  velocity  field  and 
density  is  thermal i zed.  This  continues  at  the  first  collection  radius  until 
enough  material  has  flowed  into  the  implied  central  cavity  to  erase  the 
density  gradient  to  r^.  At  this  point  a  new  r^  >  r^.  is  chosen  at  the  next 
grid  location  outside  the  previous  radius;  this  models  the  propagation  of  the 
deceleration/turbulence  front.  The  algorithm  is  continued  until  the  densest 
region  of  fluid  intercepts  the  collection  radius  and  begins  a  compression  of 
the  (now  uniformly  filled)  central  region. 

Once  the  particles  are  "collected",  the  third  phase  one  must  treat  is 
the  motion  of  the  assembled  plasma.  This  is  modeled  by  again  retaining  a 
homogeneous  7*V  and  calculating  the  dynamics  of  the  density  maximum  r^j  in  a 
manner  similar  to  that  used  in  the  run-in  phase.  In  particular,  the  x  Bg 
stress  developed  at  the  outer  boundary  is  transmitted  to  the  density  maximum 
r^^^  and  must  compete  with  the  local  internal  pressure  in  accelerating  this 
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domain  of  fluid.  The  V*^  homogeneity  then  determines  the  accelerations  and 
velocities  of  all  other  fluid  elements,  viz. 


»  °ela 
m  c  nm 


(79a) 


a  =  r_ 


(79b) 


3.  Modifications  to  the  Core-Corona  Energy  Exchanges, 
(i )  Continuous  Core/Corona  Transition. 


The  smooth  density  fall  off  with  radius  in  the  outer  regions  of  the  self¬ 
similar  core  profile  now  allows  us  to  define  the  Core-Corona  Boundary  as  the 
point  where  instability  to  ion  sound  turbulence  sets  in.  The  transition  from 
core  temperature  T^  to  coronal  temperature  T2  must  also  occur  on  some  scale 
length  \j  deterined  by  the  characteristics  of  those  modes  most  heavily  ex¬ 
cited  in  the  coronal  turbulence.  Since  ion  sound  waves  require  a  large  temperature- 
gap  T2  >  Tj  in  order  to  grow,  the  simplest  self-consistent  picture  is  one 
of  linear  temperature  transition  over  a  distance  corresponding  to  a  few 
growth  lengths  for  the  most  unstable  mode.  In  the  case  of  the  ion  sound 
waves,  this  is  a  length  on  the  order  of  several  Debye  lengths.  A  linear 

temperature  connection  and  conservation  of  current  density  j_  ,  provide  a 

c, 

simple  density  transition  n(a+XT)  =  —  n,.  Outside  the  radius  (Xy  +  a), 

I  C2  ®  J 

the  previous  isothermal  coronal  quasi-equilibrium  at  T2  is  appropriate. 

Once  this  density  and  temperature  structure  is  adopted,  it  is  possible  to 
calculate  the  heat  flux  Qj^2  particle  interchange  Fj^2*  using  the  simple 
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Krook  model  discussed  previously.  The  result  is  relatively  simple, 


Qj2  ^ 


F^.2  =  -n-a 


\2  m. 


(80) 


'  * 


’  *  <“a"ewl' 


I  "  ) 

(-  h  ^ee 

''t^ew  "I  ^ 

v\) 

'  *  <Vew>^  1 

(81) 


where  the  super-scripted  "+"  on  T,;  denotes  gradients  in  the  transition  region, 

a 

a  <r<a+\y,  T  is  %(T^  +  T2),  n  is  ^Cn-i  +  n^) ,  and  is  a  momentum  exchange 
time  for  electrons  with  unstable  waves.  The  calculation  of  F^2  ''s  an  important 
addition  to  the  coronal  temperature  derivative  calculation,  since 


^"c,e■^(W2l^'^E^-^,2*0l2  (82) 

thus, 

dt  ■  ^r,2  ^  *^12  ■  ^2^2  ’  ^^12^^!^  ’ 

with  F^2  net  electron  exchange,  weighting  by  1/Z accompanying  the  diffusion 
of  ions  across  the  interface.  These  relations  are  currently  implemented  in 
SPLAT,  and  result  in  lower  (more  physical)  coronal  temperature. 
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(ii)  Core  Heating. 


The  spatially  integrated  forms  of  the  electron  and  ion  heat  Equations 
57  and  58  are  still  useful  for  an  approximately  self-similar  core.  Since 
ohmic  heating  is  relatively  insignificant,  and  since  viscous  heating  and  pdV 
works  are  spatially  homogeneous  or  nearly  so,  and  since  (for  a  relatively 
opaque  plasma)  the  radiative  loss  is  primarily  a  surface  term,  one  can  simply 
adopt  these  previous  heating  expressions  with  only  minor  changes.  This  pre¬ 
serves  the  isothermal  nature  of  the  self-similar  aporoximation  to  the  fluid 
motion,  yet  it  allows  a  more  general  calculation  of  radiation  losses  using 
the  explicit  density  profile. 

3.  Improvements  in  Radiation  Modeling. 

At  present,  the  radiative  losses  are  calculated  on  the  basis  of  several 
new  analytic  fits  to  recent  sophisticated  steady-state  radiation  computations. 

A  frequency- integrated  volume  emission  rate  is  computed  for  the  average  plasma 
ion  density  and  electron  temperature,  to  include  bremsstrahl ung^ recombination, 
and  line  radiation  from  A1  ( K,  L)  shells.  The  line  radiation  emission  per 
unit  volume  is  fit  with  good  precision  to  the  results  of  Duston  and  Davis  (Ref.  19 
by  using  a  series  expression  of  the  form: 

n 

(S4) 

1  k 

in  (Tg/1. 200)1  .  for  T^  ^  300  ev, 

which  is  scaled  by  (fij/io^^)^,  and  subjected  to  an  attenuation  calculated  on 
the  basis  of  simple  probabil ity-of-escape  ideas  (Ref  7).  In  particular,  the 
optical  depth  is  estimated  for  a  Doppler-broadened  0.5  keV  line  as 


(Te/ 0.300)1 


for  Te<  300  ev 


74 


(85) 


T«  2.11735-17  2a{t)Rj(t); 

then  probabilities  of  escape  and  quenching 

II _ _  T  <  3 

1  +  1.86T  +  0.817t^  P  - 

n,(r)  +  1.0E21  ’ 

0.286  T  >  3 

T  (iln(1.95T))'^ 

are  calculated.  The  power  radiated  in  K  and  L  lines  is  added  to  the 
bremsstrahlung  (8r)  and  recombination  loss  (FB)  rates,  suitably  attenuated  by 
simple  optical  depth  calculations,  and  a  total  effective  emissivity 


€j(nj(r),  Tg,  Tp  r) 


1  +  Pq(Pe^  -  1) 


+  6 


Sr 


is  developed  for  any  radial  location  in  the  model  profile.  This  6-  includes 
an  additional  line  radiation  attenuation  due  to  colli  si onal  quenching  where 
the  ion  density  exceeds  10^^  cm  (Apruzese,  Ref.  8). 


With  the  net  loss  from  each  radial  location  calculated,  6.j.( r)  is 
numerically  integrated  over  the  model  profile.  The  total  radiative  loss  thus 
developed  is  then  bounded  above  by  the  equivalent  radiation  from  a  cylindrical 
black  body  at  the  current  electron  temperature. 


C.  Time-Dependent  Self-Similar  Code  SPLAT. 

1 .  Operation  and  Structure  of  SPLAT. 

The  implosion  code  SPLAT  can  be  decomposed  into  nine  distinct  sections 
covering  two  major  calculations.  It  is  equipped  with  8  load  time  parameters 
(enabling  semi-interactive  execution  and  study  of  results)  and  20  input  variables 
modeling  various  aspects  of  machine  geometry,  electrical  characteristics,  and 
initial  plasma  conditions.  In  the  present  configuration,  it  employs  a  10  radial 


fluid  grid  which  is  moved  in  a  Lagrangian  sense  and  enables  one  to  track  the 

local  neutral  plasma  position,  velocity,  and  acceleration  fields.  The  ion 
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density  is  the  fundamental  fluid  quantity  which  is  evolved,  between  the  radii 
a(t),  b(t),  in  a  self-similar  fashion.  The  ion  fluid  is  heated  through 
all  appropriate  physical  processes  and  cooled  through  expansion  and  (indirectly) 
by  radiation.  At  the  present  stage  of  development,  the  electron  density  is 
determined  from  the  homogeneous  electron  temperature  through  a  simple  approxi¬ 
mation  Z  (T  Z  A)  and  quasi -neutrality,  n^  =  Z  n..  The  electron  heating  is 
by  compression,  viscous  and  ohmic  heating,  as  well  as  by  core/coronal  heat  ex¬ 
change  and  radiative  losses. 

The  transition  from  classical  to  turbulent  current  transport  occurs  at  the 
core/corona  boundary  a(t),  and  outside  a(t)  the  appropriate  quasi-static 
coronal  pressure  balance  determines  the  electron  (and  ion)  density  profiles. 

The  coronal  electrons  are  sustained  at  T2>  T^  by  enhanced  J*E  heating  as  the 
turbulent  region  absorbs  a  larger  fraction  of  the  total  current,  while  the 
coronal  ions  are  assumed  nearly  isothermal  with  those  in  the  core  due  to  re¬ 
duced  densities  and  longer  electron-ion  relaxation  times.  Each  advance  of  the 
core  plasma  and  circuit  variables  is  initiated  with  a  calculation  of  the  inter¬ 
nal  energy  charges  in  the  coronal  region,  since  this  is  a  central  and  necessary 
ingredient  in  evolving  npr(  t) ),  T^ ,  Tj,  I .  In  fact,  all  core  plasma  changes 
T^,Tj,V(r),  and  overall  circuit  response  I  can  be  calculated  from  Eq.  50 

only  when  the  coronal  state  is  calculated.  The  sequence  of  the  time  derivative 
calculations  at  each  timestep  is  therefore  immaterial  except  for  the  constraints 
that  (i)  the  corona  variables  a  a  g  g  n  5  S  are  computed  prior  to  the  core 
response  T.| ,  Tj,  DV/Dt,I,  and  (ii)  the  radiative  losses  are  computed  prior  to 
the  core  internal  energy  changes,  allowing  a  properly  self-consistent  coupling 
of  radiative  losses  and  plasma  dynamics. 
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The  two  major  calculation  phases  of  any  SPLAT  execution  are  (i)  input 
and  normalization  (described  in  I  below)  and  (ii)  evolution  forward  from  the 
specified  initial  condition  (described  in  II-IX  below).  The  first  phase  uti¬ 
lizes  the  subroutines  NORMIT  PARCYL  PROFILE  GRID  and  SETRPC  to  initialize  a 
self-similar  ion  density  profile,  lONDEN,  on  the  domain  RGRID  with  the  charac¬ 
teristic  parameters  of  the  profile  ,  G.| ,  G^,  c^,  C^,  r^, 

stored  in  the  array  PROPAR.  The  logarithmic  gradient  3^Zn  nj(r),  (conserved 
in  strictly  self-similar  regions)  is  also  evaluated  and  stored  in  GRLNNI. 

After  a  report  on  the  initialization,  SPLAT  begins  the  second  phase  and 
evolves  this  ion  density  profile  on  its  mesh,  RGRID,  using  the  appropriate 
self-similarity  law.  Every  five  time  steps  a  short  report  is  generated  and 
every  thirty  time  steps  a  longer,  more  detailed  report  is  produced.  At  the 
time  of  assembly  on  axis,  a  simple  collection  scheme  is  employed  to  fill  the 
central  cavity  of  the  annular  self-similar  solution  and  smoothly  transform  it 
to  a  simole  cylindrical  self-similar  solution.  This  process  is  controlled  by 
the  logical  assignment  SHIFT  =.TRUE.  After  the  assembly,  the  evolution  can 
be  followed  for  as  long  as  desired,  determined  by  initial  assignments  to  the 
load  time  parameters  LIMTS  (limiting  time-step  .<  1000)  and  MAXTIM  (max.  time  £ 
250  ns) 


i.  Input  and  Normalization:  Three  namelist  files  must  be  supplied: 
GEOM  -  (a  b  i  r^)  , 

ELEC  •  (»„  Y  Lj,  2g) 

'’'-*5  =  (*  "  ■'l  "i  "  "  '> 

at  load  time,  the  execution  command  specifies 

LIMTS,  MAXTIM,  LTSP,  KCTR,  KOPF,  KVGS,  ITR,  MESH 
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as  integer  control  variables.  The  more  important  of  these  parameters  are 
defined  below: 

a,  b  outer/inner  radii  of  initial  core  self-similar  profile 


load  length 

wall  radius  /  "return  current"  radius 

estimated  lower  limits  to  initial  generator  voltage/ 
rise  time.  V  (t)  will  determine  these  in  normal 
operation.  ^ 

initial  current  through  the  load  and  voltage  across  it. 
diode  inductance  and  nominal  generator  impedance 
atomic  species  specification 

initial  mass  density  when  uni formly  distributed  on  the 
annulus  [a,b] 

initial  temperature  in  keV  for  core  electrons 
inner  radius  of  core  current  layer 

parameters  implying  the  structural  details  of  any  self¬ 
similar  profile.  4  determines  b  initially  through  Eq- 
79.  (a,  B)  provide  control  of  the  slope  of  n^Cr)  at  the 

inner  and  outer  boundaries,  n  .  sets  a  floor  on 

a, mi  n 

the  boundary  density  at  a  in  order  to  avoid  an  impossi¬ 
ble  normalization. 


LTSP  sets  convergence  criterion  for  the  time  step,  e.g., 

1000-^107  changes,  at  a  maximum,  per  timestep. 


KCTR,  KOPF  correction  factors  to  coronal  transition  layer  and 
opacity  calculations 

KVGS  control  factor  on  generator  voltage  time  scale 

ITR  initial  corona/core  electron  temperature  ratio 

MESH  number  of  points  desired  in  the  Lagrangian  core  mesh 

RGRID,  up  to  a  maximum  of  51. 
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Upon  input  of  these  namelists  and  load  time  parameters,  SPLAT  assigns  initial 
values  for  {f-jj  Tj»  1^2’  '^■j ’  *^0 ^l’^2’  ~q’  9»c,  2,  b,  a, 

b}  and  calls  SETRPC  (which  then  calls  NORMIT),  GRID,  and  PROFIL.  When 
normalization  is  successful  (the  program  stops  otherwise),  a  report  of  the 
initial  condition  is  generated.  The  general  format  of  this  report  is: 

a  r .  b 

1 

I  r 

w 

a  b  a  b 

'Jn  I  (t  =  31.  +  MAXTIM) 

g  po  max  ‘ 


n 


e 


'0 


I 


2  A  m  f 

(detailed  titles  are  omitted  here)  with  a  listing  of  RGRID,  lONDEN,  GRLNNI, 
and  PROPAR,  followed  by 
MPTC,  Sg,  g,  3 

for  the  initial  normalized  state.  (MPTC  is  the  grid  index  of  the  density  maxi- 
mum;  is  the  value  of  T(a  -b  ). 

i i .  Scalar  Time  Derivatives  Calculated. 

The  environment  of  plasma  variables  at  each  timestep  is  used  to  produce  P 

(Z2T2)  (ZiT-j)  Tj  and  I  .  The  present  version  of  SPLAT  first  calls  ZT2D0T  and 

2 

then  uses  a  simple  area  weighting  n^  to  provide  a  correction  factor  to  the  pro- 

22 

babi 1 i ty-of-escape  formalism  if  nj>10  where  quenching  is  severe.  In  those 

22 

time  frames  for  which  n^  <10  for  all  radial  grid  points  this  additional  correc- 
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tion  is  not  involved  and  the  mean  density  fij  =  n^/Z,  is  used  to  evaluate 
a  volume  emission  rate  modified  by  simple  probability  of  escape  factors  for 
a  plasma  of  radius  a  (t)  and  an  optical  depth  computed  from  Tj,  n^.  Once 
is  calculated,  ZTIDOT,  TIDOT  and  IDOT  are  invoked  and  an  updated  set 
Ti ,  J2,  Tj,  I  is  available. 

i i i .  Velocity  (V),  Acceleration  (DV/Dt),  and  (D^V/Dt^)  Fields 

The  current  values  of  PROPAR,  T,  Tj,  I,  I  are  used  to  produce  the  fields 

required  to  push  the  Lagrangian  mesh  RGRID.  For  example,  the  use  of  G.| ,  G2 

defines  a  field  V  s  DR/Dt: 

V(R(t)) 

iv.  Diagnostic  Update 

If  the  calculation  is  about  to  advance  to  a  time  step  index  (TSI)  divis¬ 
ible  by  30,  then  a  vector  of  diagnostics  is  prepared  for  output.  This  diag¬ 
nostic  array  consists  of 

^^a  ^2  '^s.p.c.  '^corona^'^w  ^core,  max^’ 

the  electric  fields  on  the  core  surface  and  outer  plasma,  the  number  of  sound 

wave  transits  across  the  corona  in  the  time  for  a(t)  to  change  by  10",  the 

corona  radius  relative  to  the  wall  radius,  and  the  maximum  current  admissible 

in  the  core  surface  if  this  zone  is  to  rema i n  classical. 

V.  Time  Step  Determination  ("Update  of  Pushdown  Stack") 

In  this  section  the  rates  of  change  are  surveyed  and  a  time  step  is  chosen 
on  the  basis  of  several  criteria: 

(i)  fractional  changes  in  the  most  rapidly  varying  scalar  quantity  of 
the  set  (T^ ,  Tp  T^,  I) 


+  G,(t) 


TTtV 
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(ii)  possible  interactions  of  ,  Tj  to  quench  temperature  "sloshing" 

(iii)  possible  zeros  in  T2,  Q12  to  quench  core/corona  energy  "sloshing" 

(iv)  possible  zeros  in  I  to  quench  corona/circuit  energy  "sloshing" 

(v)  velocities  and  accelerations  on  the  mesh. 

vi .  Advance  of  Vector  and  Scalar  Fields  (includes  temperatures  and  current) 

The  multi-level  integrator  advances  {T.| ,  T^,  T2,  I},  while  the  mesh 
RGRID  is  advanced  through  a  third  order  Taylor  series.  If  the  collection  pro¬ 
cess  is  active  in  the  interior  region  no  mesh  advance  takes  place  on  those 
domains. 

vii.  Update  Plasma  Density  Variables,  etc. 

Here  the  new  plasma  variables  are  transformed  to  those  intermediate  quan¬ 
tities  required  to  calculate  the  new  rates  of  change.  The  ion  density  array 
lONDEN  is  evolved  according  to  the  self-similar  evolution  rule,  and  the  modi¬ 
fication  of  3^  tn  nj  in  the  current  carrying  layer  is  calculated. 

viii.  Collection  in  Central  Cavity  when  required 

The  collection  of  matter  flowing  past  an  implicitly  time  dependent  radius, 
ERC,  is  computed  and  interior  density  profiles  and  gradients  are  smoothed.  The 
collection  persists  as  long  as  matter  can  be  drawn  onto  the  central  plasma  body 
by  a  negative  VP  force,  but  it  is  interrupted  (HOLD  =  .TRUE.)  if  the  local 
flo,v  field  is  not  directed  toward  the  central  plasma. 

ix.  Housecleaninq  and  Reporting 

Here  preparation  for  the  next  time  step  are  completed  by  a  call  to 
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SETRPC,  and  a  report  is  constructed  if  required.  The  shorter  report  output 
every  5  time  steps  has  the  format 
a  r^.  b  a  b  C2 

I  Vg  R* 

"e  ^2  ''’l  '’rad  ^12 

t  dt  TSI 

a  b  Cl  B  Qi2-W2 

RGRID 

lONDEN 

The  longer  report  presents  several  more  details  of  the  evolution,  information 
permitting  a  cross  check  of  the  precision  of  the  integration,  the  array  DIAG, 
and  details  of  the  coronal  state,  CPV. 


2.  Description  of  Subroutines  and  Functions 

(1.)  SETRPC  -  Set  radial  profile  coefficients 

•  argument  sequence:  LOGPRO,  n,,  N.,  Z, ,  m,  T,  fi, ,  n,,  E' ,  a,  r . , 

oil  I  j.  d  1 

•  * 

b,  a,  b,  a*  b 

•  output  to  COMMON  BLOCK:  DENSI/PROPAR/  and  (if  LOGPRO  =  .FALSE. )  the 
argument  "a"  is  overwritten  by  the  subroutine  NORMIT 

•  operational  summary:  If  the  initial  profile  has  not  been  calculated 
and  normalized,  SETRPC  prepares  the  common  block  ZNARG/A$/  for  use  by 
NORMIT  and  then  calls  NORMIT  to  adjust  "a"  until  the  profile  contains 
a  specified  linear  density  Ny.  If  the  profile  has  been  calculated  and 
normalized  previously  (LOGPRO  =  .TRUE. ) ,  then  only  the  elements  a^,  a.j , 
G.|,  G2,  (i.e.  (PROPAR  (1,  2,  3,4))  are  updated  using  the  current 
values  of  a,  b,  a,  b,  a,  B. 
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( 2.)  NORMIT  -  Normalize  It. 


•  argument  sequence:  a  r^  b  a  b  c-j  a  B 
input  common  blocks  ZnARG/A$/,  CORO/CPV/ 

•  output;  overwriting  of  argument  "a"  ;  and  to  COMMON  BLOCK;  CORO/CPV/ 

to  update  3,  g,  i.e.,CPV(5),  CPV(6)  ,  after  a  normalization  calculation. 

•  operational  summary:  NORMIT  solves  the  equation  (  77)  for  that  value 
"a"  which  implies  a  specific  total  number  of  elections  per  unit 
length  in  the  core  plasma.  AS(1)  contains  this  total  number  upon 
any  call  to  NORMIT.  The  iteration  is  continued  until  convergence  to 
four  significant  figures  in  the  value  "a"  is  achieved  or  until 
twenty  iterations  have  occurred.  Usually  three  to  seven  iterations 
are  quite  sufficient  to  isolate  the  root;  each  iteration  is  reported 
as  a  triangle  of  points  enclosing  the  root  and  a  slope  for  the  hypote- 
neuse  of  this  triangle.  The  value  of  the  denominator  on  the  right  of 
Eq.  77,  with  the  numerically  integrated  contribution  to  it,  is  reported 
along  with  the  values  generated  in  PROPAR  as  the  iteration  proceeds. 

If  convergence  does  not  occur  the  final  state  of  the  calculation  is 
reported  for  diagnosis.  Since  ng(r)  involves  the  special  function 
U,(xR,  u,  a)  NORMIT  calls  PARCYL  whenever  such  functions  are  required. 

(3.)  PROFIL  -  Profile 

•  argument  sequence:  rr],[nj(r)]  ,  [a^2,n  n^],  z^,  ii-i ,  m,  f,  array  dimensions 
input  common  block:  DENSI/PROPAR/ 
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•  output:  overwriting  of  the  array  arguments  Cnj(r)3  nj(r)] 

•  operational  summary:  PROFIL  calculates  the  explicit  self-similar 
profile  for  an  array  of  values  r  using  the  parameters  stored  in 
PROPAR.  It  calls  PARCYL  whenever  the  special  function  U-(xR5  -) 
is  required.  It  stores  these  values  in  nj(r),  the  array  lONDEN 
from  SPLAT.  It  also  computes  3^  in  n^  and  stores  these  values  in 
the  array  GRLNNI  from  SPLAT. 

(4.)  PARCYL  -  "Incomplete  Parabolic  Cylinder  Function'' 

•  argument  sequence:  [U]  [X]  R,  u,  a,  :rray  dimensions, 

•  output:  overwriting  of  the  array  argument  [U] 

•  op-?rational  summary:  PARCYL  evaluates  the  special  function 

U_(/R,  u,  a)  by  numerical  integration.  The  accuracy  of  the  integration 
is  determined  by  the  variable  DG  (dimension  of  the  grid),  computed 
externally  on  the  basis  of  -  and  R  and  bounded  above  by  DG_.^  =  301. 
Seven  figure  accuracy  is  easily  achieved  on  a  2  [lO  +  8  („  R)]  +  1 
mesh  for  x  $1 . 

( 5. )  GRID  -  Mesh  Generator 

,  argument  sequence:  Cr]  r . ,  r^.^,  I^  ,  points 

i 

•  output:  overwriting  of  array  argument  [r] 

•  operational  summary:  GRID  produces  a  sequence  of  equally  spaced  r 

values  between  r  •  and  r.,  and  inserts  four  points  between  r.  and  r  . 

The  grid  then  consists  of  (points  -  5)  values  on  rr„.  ,  r.]  and  five 

mm  1 

values  on  (r^. ,  .  for  a  total  of  (points  -  1)  zones. 
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(6.)  ZT2D0T  -  Corona  Internal  Energy  Time  Derivative 


•  argument  sequence:  n^,T^,T2,  I,  E^,  A 

a  r.  0„  Z,  T,  E'  n.  3  2n  n, 

1  p  1  T  I  a  I, a  r  I, a 

input  from  common  block  CORO/CPV/:  S  g  _ 

0  I  I  jL, 

CPV  (3,  6,  13,  10) 

•  output:  to  common  block  CORO/CPV/  T2  ^  "  n^  P^  2  Q-i2  ^12 

CPV  (1,  2,  7,  8,  9,  11,  12,  14) 


•  operational  summary:  ZT2D0T  calculates  the  coronal  equilibrium  para¬ 
meters  3,  n,  5  (and  N-  on  the  first  time  step).  The  core  (corona 

I  *L 

density  and  temperature  transition  is  smoothly  connected  by  conserving 
over  a  narrow  region  '•lj.  on  the  order  of  a  few  wave  growth 
lengths.  Using  the  explicit  density  gradient  supplied  from  the  core 
and  differencing  over  to  estimate  the  density  and  temperature 
gradients  at  the  inner  surface  of  the  corona,  ZT2D0T  distributes  the 
net  energy  gain  or  loss  over  all  Nj  j,  atoms  in  the  coronal  equilibrium 
and  generates  values  for(T2,  ct  ,  Q.|2.  ^12^  ’  Krook  model  is  applied  ! 


consistently  in  all  exchange  calculations  leading  to  Q.|2  (kinetic 

energy  transfer)  and  F^2  (particle  transfer).  The  far  coronal  field 

E  is  used  to  calculate  the  J-E  heating  and  the  functions  RAD  and 
z  z  z 

CORFAC  are  used  to  compute  the  radiative  loss  from  the  thin  corona. 
(7.)  CORFAC  -  Corona  Radiation  Factor 


•  argument  sequence  (r,  n  »  5.  term  limit) 

•  output:  to  function  nami  ^CORFAC 
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•  operational  summary:  CORFAC  evaluates  the  geometric  weighting  fac¬ 
tor  in  the  net  radiative  energy  loss  from  a  thin  isothermal  radi¬ 
ator  with  a  density  profile  n^{r)  corresponding  to  the  coronal 
equilibrium  configuration  given  at  any  time  step  by  r,  n,  o  . 

This  evaluation  is  accomplished  either  by  hypergeometric  series  or 
analytic  approximation  where  appropriate. 

( 3. )  RAD  -  Radiation 

•  argument  sequence:  (n  Z  A  T  <.  7.  0  ) 

0  0  1  p 

•  output:  to  function  name,  RAD 

•  operational  summary:  RAD  calculates  a  frequency  integrated  volume 
emission  rate  for  a  plasma  soecified  by  (n^,  z,  Tj,  A).  The 
Bremsstrahlung  and  recombination  losses  are  presumed  thin,  while 
the  K,  L  shell  line  radiation  losses  are  calculated  from  an  analytic 
fit  to  the  emission  rate  results  of  Duston  and  Davis  (the  function 
DAVSR)  and  then  attenuated  using  the  probabil ity-of-escape  opacity 
model  due  to  Apruzese.  The  arguments  n^,  z,  I,  Tj,  Op  control  this 
opacity  calculation;  I  plays  the  role  of  the  physical  depth  or  dimen¬ 
sion  of  the  radiating  plasma,  i=0  results  in  a  thin  calculation.  Argu¬ 
ments  n^/Z,  Tj  are  used  to  calculate  an  appropriate  optical  depth  and 

Op  allows  any  external  calculation  of  opacity  corrections  to  be  included 
in  the  attenuation. 

(9.)  DAVSR  -  Line  Emission  Due  to  Duston  and  Davis 

•  argument  sequence:  n^  T^ 

•  output:  to  function  name, DAVSR 
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(9. )  (cont'd) 

•  operational  summary:  This  function  fits  the  K,  L  shell  A1  line 


emission  rates  calculated  by  Duston  and  Davis  over  the  temperatures 

1  Q 

range  =  (0.010,  9.000)  kev  and  scales  this  emission  by  (nj/l.OXlO  ) 
The  fit  is  good  to  within  Q%  over  the  entire  range  of  temoeratures . 

( 1 0 . )  ZTIDOT  -  Core  Internal  Inerqy  Time  Derivative 

•  argument  sequence:  n^  T^  Tj  r  *^1 

a  b  c^  a  N,  2^  A  a  b 

•  output:  to  function  name, ZTl DOT 

•  operational  summary:  ZTIDOT  calculates  the  core  heating  due  to  com¬ 
pression.  [p(T*V.),  ohmic  dissipation  ,  viscosity 

core/corona  energy  exchange  Q.|2.  total  radiative  energy  loss  P^,  and 
ion-electron  energy  exchange.  The  net  energy  gain  or  loss  is  distri¬ 
buted  over  the  conserved  number,  of  atoms  in  the  core  plasma.  All 
radiative  loss  rates  are  computed  externally  using  RAD  and  DAVSR  and 
allowed  to  enter  self-consistently  into  the  internal  energy  evolution 
of  the  core  plasma.  Under  conditions  of  particularly  violent  compres¬ 
sion,  a  heating  report  is  printed  which  shows  the  relative  importance 
of  the  various  heat  sources  during  the  time  of  interest. 

(11.)  TIDOT  -  Core  Ion  Heating 

•  argument  sequence:  rig  T.j  Tj  Z.j  C.|  Ny  A 

a  b  a  a  b 

•  output:  to  function  name, TIDOT 
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(11.)  (cont'd) 


•  operational  summary:  TIOOT  calculates  the  core  ion  heating  due  to 
compression  [p(7>^)] ,  electron-ion  heat  exchange,  and,  if  desired, 
an  enhanced  viscosity  Cn*  under  conditions  of  rapid  compression 

(a|>  1.5c-].  The  viscosity  enhancement  is  smoothly  connected  to  the 
classical  value  as  !aj  falls  below  the  current  sound  speed  c-j .  Under 
conditions  of  particularly  violent  compression,  the  heating  report 
started  by  ZTIDOT  is  completed  by  TIDOT. 

(12.)  IDOT  -  I,  Calculation  from  the  Circuit  Relation 

•  argument  sequence:  a  a  i  r^  I  R*  Lq 

adgg  Sn  6  V 

b 

•  output:  to  function  name,  IDOT 

•  operational  summary:  This  function  constrains  the  time  derivative 

I  to  be  given  by  the  complete  circuit  equation  (50)  with  all  coronal 
contributions  explicitly  included.  The  driving  voltage  is  supplied 
by  Vg  and  all  other  input  is  derived  from  the  results  o^  ZTZDOT,  and 
preparatory  calculations  in  SPLAT  from  current  values  at  a  given  time 
step.  When  a  decrease  in  current  is  imminent,  a  report  of  IDOT  DATA 
is  generated  to  allow  diagnosis  of  the  factors  controlling  the  "current 
notch" . 


(  13.)  SET5T  -  Set  Time-steF 


•  argument  sequence:  a,  a,  a,  a’,  b,  c-j ,  convergence  and  history  information 

•  output:  to  function  name, SETST 
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•  operational  summary:  SETS!  surveys  the  most  rapidly  varying  elements 
of  the  evolution  and  chooses  the  time  step  required  for  a  predeter¬ 
mined  level  of  accuracy  (the  LISP  in  SPLAT)  in  the  most  rapidly 
varying  of  (a,T^,  Tj High  frequency  exchanges  of  energy  are 
anticipated  and  quenched  to  allow  the  time  step  to  remain  above 
5.0  E-11  sec. 


(lA.)  WSPUSH  -  Weights  for  the  Variable  Time  Step/Mul ti -time  Level  Integrator 

•  argument  sequence:  C^JCh] 

•  output:  overwriting  of  argument  array  [W] 

•  operational  summary:  WSPUSH  computes  the  weight  factors  for  advancing 

Ti  .Tj.T^’^  through  the  leap-frog  integrator  =  ( 2.^  ^i  2 

The  variable  mesh  of  previous  time  steps  QhJ  is  required  in  calcu¬ 
lating  W^. 

(15.)  ETA  -  Classical  Resistivity, 

•  argument  sequence:  (n^  T^  Z) 

•  output:  to  function  namejETA 

•  operational  summary:  ETA  evaluates  the  resistivity  n  given  by 

n  =  a  Z  ^n  >  T‘^^‘ 

XX  e 

(16.)  V  -  Generator  Voltage  Model 

•  argument  T  in  nsec 

•  output:  to  function  name^V 

•  operational  summary:  a  piecewise  quadratic  fit  to  Vg(t)  due  to  J. 
Katzenstein. 
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3.  SPLAT  Program  Flow 


READ  3  namelists  GEOM,  ELEC,  PLAS,  which 
jSet  generator  parameters,  initial  plasma  state. 

I^nitialize  a,  a,  a,  b,  b,  B;  ni,T.j,  Ip  3,  N. 

Normalize  an  initial  orofile  to  enveloD  a  soecified  mass;  'Jm,, 

A  1 


Call  subroutine  SETR(adial )P(rofile)C(oefficients) , 

which  call  subroutines  NORMIT  and  PARCYL  to  adjust  profile 
width,  and  evaluate  Gp  2.,,  'Up  (e.g.  31),  the  elements 
(e.g.  50)  output  through  the  common  block  PROPAR  (Profile 
PARameters ) . 

Initialize  the  arrays  lONDEN,  R6RID,  and  GRL.'iNI  (gradient  of  loga¬ 
rithm  of  ion  density)  usino  the  subroutines  GRID,  PRCFILE,  PARCYL 

] 

Report  the  initial  state 


Calculate  corona  oarameters  using  subroutine  ZT2D0T,  and  outout  through 

the  common  block  C(orona)  P(arameter)  V(ector)  the  variables  (t2,r,  i, 

Qp,  P^  2’  2T200T  calls  RAO  (which  calls  DAVSR)  and  function 

CORFAC  to  effect  a  thin  calculation  of  corona  losses. 
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Calculate  the  relative  opaque  volume  of  core  plasma,  cf.  Eq.  (86),  for 
use  with  probabi 1 ity-of-escape  loss  calculations, 

[Calculate  core  radiative  loss  using  RAD  and  DAVSR . 

ZTIDOT  calculates  ,  using  the  present  heat  flow  radiated  power 
and  fluid  plasma  parameters. 


I 


TIDOT  calculates  T,  for  the  present  fluid  plasma  parameters. 

i  . 

Calculate  r^  from  the  present  flow  field  in  common  block  PROPAR,  and 
the  appropriate  g. 

Calculate  the  total  I  using  the  complete  core/corona  circuit  equation  and  the 
present  values  of  g,  g,  a,  a,  3,  5,  n,  R*,  using  the  function  IDOT  . 

3 

Prepare  to  advance  the  Lagragian  mesh  RGRID  using  the  sel f-s'.mi  1  ar  velocity 
fields  AVSSC  (-,  3),  and  acceleration  fields  AVSSC  (-,  4). 

AVSSC  (-,  5)  calculated  from  PROPAR. 


3 


[Prepare  the  array  OIAG  if  required  (every  30  time  steps';  tor  printout 


Push  down  the  multi-time-level  stack  of  t.j  tj  I  T.|  "'"t  "^2  ^  , 
the  sequence  of  heat  flows  Q12S,  and  the  sequence  of  time  steps  HS 

3 

Inventory  the  rates  t^/T^,  tj/Tj,  t^/T^,  l/l  and  calculate  the  t 
step  using  SETST(ep)  for  DT. 

3 


ime 


Advance  the  portion  of  RGRID  not  involved  in  collection  and  predict  new 
velocities  AVSSC(-,2)  and  accelerations  AVSSC(-,1)  for  the  zone  boundaries 
a  and  b. 
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Elall  W$PUSH  to  calculate  the  weights  for  the  multi -time-level 
integrator,  and  advance  Tj,  T2,  I- 

N  — the  integrator  has  produced  variations  beyond  tolerance 
Y 

[Restore  and  use  Euler  method 

1 

——^Advance  corona  particle  number 

Advance  the  clock  (T)  and  T(ime)  S(tep)  I(ndex) 


Update  a,  b,  a,  b,  a,  from  the  mesh 
Calculate  ,  Z2,  ri  ^ R*,  £ 

Update  lONDEN  using  the  s/s  evolution  rate. 

Update  the  array  GRLNNI  in  the  current-carrying  layer  (which  is  not 
strictly  self-similar)  using  differencing  estimates. 

'J - |lf  collection  of  the  central  plasma  is  needed  or  ongoing  (SHIFT  =  TRUE.) 

V 

Collect  ions  flowing  past  a  small  specified  raoius  ERC  until  they 
uniformly  fill  the  central  cavity 

Maintain  strict  s/s  soTn  outside  ERC  using  current  boundary  values 
of  r^,  r^ 

— — -^Update  PROPAR 


Report  plasma  state  as  appropriate 
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MICROCOPY  resolution  TEST  CH4RT 


Check  for  ending  conditions - MAXTIM,  LIMTS,  refrigerati ve 

collapse. 


Y 

[STOP 
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4.  Advantages  and  Disadvantages 

Advantages: 

(i)  Economy  of  computation  (Ic/time  step) 

(ii)  Relative  simplicity  in  treating  new  physics  (zone/module 
structure ) 

(iii)  No  arbitrary  multipliers  on  transport  coefficients 

(iv)  A  class  of  initial  conditions  of  specific  form  that  are 
subject  to  experimental  constraint,  reducing  the  dimension 
of  the  free  parameter  space 

(v)  Potential  generalization  to  2  and  3  dimensions  with  even 
greater  relative  computational  economy. 

Pi sadvantaqes : 

(i)  Breakdown  of  self-similar  motions  upon  collapse  should  be 
treated  with  more  general  fluid/kinetic  theory,  but  the  self¬ 
similar  regions  allow  simple  dynamic  calculations  and  free 
the  more  delicate  methods  to  concentrate  on  the  special  cases. 

This  implies  less  computation  for  any  fixed  precision  or  more 
precision  for  a  fixed  level  of  computational  effort. 

(ii)  The  very  stiff  set  of  non-linear  differential  eguations  requires 
delicate  integration  methods.  Particularly  if  extended  into 
domains  of  sharp  gradients  and  many  orders  of  magnitude  in  den¬ 
sity,  the  tracking  of  widely  separate  rates  for  various  energy/ 
momentum  exchanges  is  the  most  difficult  aspect  of  the  calculation. 
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V.  SCALING  AND  GUIDELINES  BASED  ON  CORE-CORONA  MODEL  AND  CODES 

A.  Interpretation  of  Code  Results. 

1 .  Initial  Conditions. 

The  self-similar  initial-condition  relations  discussed  above  allow 
a  wide  class  of  initial  core  states,  and  there  is  little  experimental  infor¬ 
mation  to  provide  more  than  a  rudimentary  guide  to  their  selection.  The 
parameters  of  the  self-similar  profile  a^,  G^,  G2,  m,  T,  and  E^  are 
chosen  for  the  code  runs  on  the  basis  of  reasonable  criteria,  but  are  not 
as  yet  intended  to  duplicate  any  particular  laboratory  implosion. 

In  particular,  one  wishes  an  initial  profile  exhibiting:  (i)  confine¬ 
ment  to  a  narrow  annulus,  (ii)  expansion  velocities  on  the  order  of  the 
sound  speed,  (iii)  expansion  velocities  compatible  with  initial  temperatures 
(of  order  lOeV),  (iv)  initial  currents  compatible  with  the  generator  voltage 
and  network  parameters,  and  (v)  a  specific  total  mass.  The  term  "compatible" 
means  that  in  a  practical  sense  one  must  search  the  parameter  space  a  bit 
with  each  fresh  profile  constructed,  usually  adjusting  a,  b,  I  in  order  to 
achieve  satisfactory  confinement.  The  choice  of  r^. ,  the  extent  of  radial 
penetration  of  the  driving  electric  field  prior  to  freezing  these  fields  to 
the  fluid  flow,  is  also  relatively  unguided  by  direct  experimental  input. 

At  the  temperatures  considered,  the  resistivity  is  so  low  as  to  force  any 
statement  of  to  be  a  statement  about  the  early  time  development  of 

melting,  vaporization,  expansion  and  ionization.  None  of  these  transitions 
are  intended  to  be  a  part  of  this  model  and  thus  their  end  results  are 

treated  as  simple  input  data.  For  example  at  32  ev,  a  skin  depth  of  8.4  x 

.3 

10  cm  is  well  within  reason  for  times  as  late  as  45  ns  into  the  implosion 
phase  where  the  self-similar  model  could  reasonably  be  expected  to  accommodate 
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t  the  physics.  The  selection  of  the  present  initial  profiles  as  physically 

appropriate  rests  on  the  suppositions  that:  (i)  the  early  fluid  development 
passes  through  a  state  of  large  viscosity  (relative  to  the  classical  plasma 
value)  and  (ii)  such  viscous  stresses  act  to  bring  the  density  and  flow- 
field  profiles  into  the  forms  used  here.  Whatever  the  detailed  time  develop¬ 
ment  in  the  early  phases,  the  restriction  to  such  a  class  of  profiles  is  a 
relatively  weak  one.  It  implies  generation  processes  which  provide  a  con¬ 
tinuum  of  initial  values,  e.g.,  r^,  a,  b,  and  which  admit  further  experi¬ 
mental  investigation. 

The  fixed  initial  mass  is  taken  to  be  that  of  a  24-wire  array  of  I. 5 
mil  Al  wires  3  cm  long.  This  is  used  by  the  normalizing  algorithm  discussed 
above  and  is  a  conserved  quantity  in  the  calculation.  The  structure  of  two 
such  self-similar  profiles  is  shown  in  Fig.  (3a, b)  for  initial  temperatures 
of  1 5  eV.  The  smaller-radius  profile  initially  confines  most  of  the  mass 
between  0.60  cm  and  0.76  cm,  the  larger-radius  one,  between  0.92  cm  and  I. 09 
cm.  Both  profiles  contain  I. 5023  x  lO^^  ions  per  cm  of  discharge  length. 

The  initial  penetration  of  the  current  layer  is  indicated  at  0.60  cm  and  0.95 
cm,  respectively. 

2.  Classes  of  Trajectories. 

In  the  studies  detailed  below,  one  can  identify  three  general 
classes  of  implosion  trajectories,  as  defined  by  the  position  of  the  density 
maximum  r^(t)  in  the  self-similar  profile.  The  first  is  a  simple  refrigera- 
tive  collapse  (SC),  in  which  the  implosion  picks  up  energy  from  the  generator, 
builds  inadequate  retarding  pressure,  and  quickly  radiates  this  energy  away 
as  it  is  thermalized.  Under  these  conditions,  the  plasma  is  smoothly 
driven  into  the  origin  by  J_xB_  forces  without  evidencing  any  retarding 
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Figure  3  a,b.  Two  "initial"  ion  density  profiles  for  differti  fc  SPLAT-ID 
code  runs.  Initial  temperatures  were  taken  to  in,  15  eV. 

The  assumed  initial  current-penetration  radiu-  ,  r. .  is  shown 
for  both.  ’ 


pressures  large  enough  to  cause  deceleration.  The  second  type  of  trajectory 
pauses  and  collapses  (PC)»  where  the  rapidity  of  heating  allows  retarding 
pressures  to  build  up  and  slow  the  implosion  velocity  by  a  factor  of  2  or 
more.  In  this  class  the  radiation  output  becomes  slightly  broader  in  time 
and  larger  temperatures  are  seen  (nearly  coincident)  with  the  radiation  peak. 
A  third  type  of  trajectory  involves  more  significant  bounce  before  collapse 
(BC),  and  runs  with  higher  peak  temperatures  and  faster  generator  risetime 
fall  into  this  latter  class.  The  bouncing  occurs  in  the  zero-0  model  also, 
and  is  more  asymmetric  in  time  if  the  outward  pressure  gradients  after 
implosion  are  more  pronounced;  the  bounce  is  held  to  smaller  radii  only  if 
the  energy  feeding  expansion  is  removed  rapidly  enough  by  the  radiative 
losses.  Bounce  trajectories  without  collapse,  although  not  seen  in  the 
present  study,  can  be  generated  for  weak  or  rapid  generator  pulses. 

The  density  and  temperature  trajectories,  which  determine  the  radiation 
pulse  width  and  amplitude,  are  similar  for  the  SC  and  PC  paths,  and  distinct 
from  those  of  the  BC  paths.  In  the  case  of  the  former,  all  densities  on 
the  profile  increase  monotonically  with  time  after  the  Initial  expansion 
phase  of  the  annulus.  The  temperatures  peak  just  prior  to  peak  assembly 
on  axis,  shortly  before  the  maximum  radiative  losses.  In  the  BC  case,  the 
densities  rise  until  the  bounce  occurs,  after  which  there  is  a  smooth 
decay  and  subsequent  re comp res si  on.  The  temperatures  are  observed  to  rise 
and  fall  smoothly  through  the  initial  bounce,  only  to  rise  again  from  small 
values  C-- 30  eV)  with  the  final  compression  and  collapse. 

In  all  cases,  some  time  asymnetry  of  the  radiation  pulse  occurs  because 
the  rise  time  is  determined  by  the  plasma  compression  time  scales,  but  the 
decay  is  determined  mainly  by  the  large  radiative  loss  rate  which  is  very 
sensitive  to  even  the  slowed  compression  or  expansion  of  the  pinch.  The 
mix  of  line  (K,  L  shell)  and  continuum  radiation  tends  to  favor  the 


continuum  at  the  peak  of  emission  because  the  high  densities  n^-^  3  x  10^’  imply 
very  high  line  opacity.  In  the  absence  of  motional  broadening,  a  mix 
favoring  line  emission  will  be  obtained  for  milder  compression,  i.e., 
bounce  or  pause  at  earlier  times  and  larger  radii,  or  for  wanner  plasma. 

Larger  initial  temperatures  in  the  annulus,  perhaps  50  ev  -  70  ev,  tend  to 
aleviate  this  behavior,  but  many  alternative  factors  can  play  a  significant 
role. 

The  time-integrated  yield  (integrated  over  all  photon  energies)  is 
sensitive  to  the  maximum  implosion  speed  achieved,  whether  by  faster  rise 
times  in  the  generator  (larger  forces),  or  by  a  longer  distance  over  which 
to  act  (greater  initial  radii).  In  the  case  of  PC  trajectories  there  will 
tend  to  be  a  short,  cold  after-pulse  following  the  main  radiation  pulse, 
caused  by  the  "last  gasp"  of  the  system  as  collapse  occurs.  The  contribu¬ 
tion  to  the  yield  from  a  single  such  pulse  is  small,  but  if  the  initial 
bounce  point  is  approached  slowly  there  may  exist  implosion  trajectories 
which  hold  at  moderate  radii  for  some  time,  radiating  more  slowly  and  bounc¬ 
ing  gently,  perhaps  many  times.  Under  these  conditions  the  collapse  could 
be  stretched  out  and  the  contribution  to  the  yield  due  to  this  late  time 
behavior  would  increase.  However,  this  contribution  is  at  photon  energies 

considerably  less  than  1  keV. 

3.  Illustrations  of  the  Dynamics. 

The  typical  SC  trajectory  is  shown  in  Fig.  4a,  starting  from  the 
initial  state  of  Fig.  3a.  The  largest  implosion  speed  of  the  density  maximum 
is  about  1  x  10’  cm/sec  and  the  peak  temperature  on  axis  is  338  eV. 

The  current  notch  due  to  motional  impedence  is  not  pronounced. 

The  typical  PC  trajectory  is  shown  in  Fig.  4b,  starting  from  the  initial 
state  of  Fig.  3b.  The  largest  implosion  speed  of  the  density  maximum  is  about 
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2  X  10^  cm/sec  and  the  peak  temperature  on  axis  is  378  eV.  The  current 
notch  is  relatively  weak,  subtracting  about  10%  of  the  previous  maximum 
driving  the  implosion;  but  the  current  trace  recovers  to  100%  of  this 
previous  current  maximum.  Even  though  the  generator  voltage  continues  to 
rise,  the  "impedance"  of  the  array  increases  significantly  as  the  radius 
decreases.  In  the  first  of  the  cases  illustrated,  the  core  current  frac¬ 
tion,  8.  is  approximately  0.7,  but  the  coronal  temperature  never  exceeds 
10  keV.  In  the  second,  8  is  bounded  above  by  0.9,  and  T2,  again  by  10 
keV.  Both  cases  were  run  with  good  initial  current  penetration  into  the 
core  plasma. 

Figures  5  and  6  illustrate  the  effect  of  variations  in  generator 
rise  time.  The  yield  dependence  on  the  generator  timescale  is  shown  for 
several  cases  with  the  initial  conditions  (Figures  3,a,b,  etc.)  and  clearly 
indicates  that  too  rapid  a  coupling  degrades  the  yield  at  some  point. 

The  generator  rise  times  were  varied  using  a  simple  multiplication  on  the 
time  scale  of  the  driving  function  V(t),  i.e.,  V(t)-^V(at  +  31.0  nsec), 
a  =  0.5  to  1.6  in  intervals  of  0.1.  (Figure  5.)  The  yield  degrades  for 
the  faster  time  scales,  because  the  generator  then  can  heat  the  plasma  earlier 
in  assembly  phase  and  cause  it  to  expand  before  reaching  as  high  a  density 
(Fig.  6).  This  provides  higher  temperatures,  but  generally  lower  densities 
and  keeps  the  overall  radiated  energies  smaller. 

Note  that  the  radiated  energy  plotted  in  Fig.  6  (and  the  radiated 
power  in  Figs.  4a  and  4b)  are  frequency- integrated  values,  without  any 
cutoff,  and  do  not  represent  experimentally  measured  radiation  passed  first 
through  filters.  The  radiation  package  of  the  code  is  presently  being  recon¬ 
figured  to  ^  ,  ot  uc  separately  in  two  frequency  "bins". 
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gure  5.  Generator  voltage 
pulse) . 


Frequency-integrated  radiative  yield  vs  generator  rise-rate 
multiplier  u,  for  three  initial  configurations; 


4.  General  Comparison  of  Zero-D  and  1-D  Formulations. 


The  use  of  a  continuous  profile  for  the  core  density  provides  an 
important  correction  to  the  dynamics  of  the  implosion.  It  is  clear  a  priori 
that  the  mechanism  for  coupling  generator  fields  to  radial  fluid  momentum 
is  a  major  determinant  of  the  dynamics,  but  it  is  obvious  after  some 
experience  that  this  is  the  single  most  important  facet  of  the  problem. 

The  use  of  explicit  pressure  gradients  and  the  propagation  of  external 
JxB  stresses  to  the  local  density  maximum  allows  the  model  plasma  to  respond 
to  the  "magnetic  piston"  by  a  smooth  inward  drift  reflecting  the  appropriate 
inertia  and  by  relative  compressions/expansions  reflecting  the  density  pro¬ 
file.  The  simple  spatially  integrated  model  only  allows  the  average  inertia 
to  enter  and  accelerates  the  inner  surface  independently  of  the  pressure 
gradient.  This  results  in  larger  implosion  speeds,  higher  temperatures, 
and  shorter  radiation  pulses  as  can  be  seen  from  the  dynamics  presented  in 
previous  sections.  On  the  other  hand,  the  smooth  ID  version  prevents  such 
large  transfers  of  radial  momentum. 


The  second  most  important  correction  is  the  opportunity  to  do  a 
detailed  calculation  of  the  radiative  losses.  From  the  simple  picture,  we 
concluded  that  proper  opacity  calculations  are  very  important  in  the  model 
dynamics  because  energy  not  radiated  remains  to  provide  more  retarding  pressure 
(directly  or  indirectly}  and  slows  the  time  scale  of  tne  radiation  pulse. 
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The  third  correction  is  embodied  in  the  core/corona  heat  exchange 
calculation  discussed  above.  This  energy  exchange  plays  a  major  role  in 
determining  B  and,  thus,  in  modifying  the  ^xB^  force  driving  the  core  implo¬ 
sion.  In  contrast  to  previous  calculations,  the  present  ID  version  allows 
for  particle  interchange  at  the  core/corona  boundary  and  provides  a  cooler 
corona,  less  susceptible  to  "exploding".  The  cooler  corona  now  tracks  the 
core  temperature  more  closely  and  is  less  extended  spatially  than  in  the 
previous  version  of  this  model.  In  fact  the  corona  now  resembles  a  thin 
zone  of  turbulent  current  transport  capable  of  accepting  essentially  all  of 
the  current  required  by  the  generator  field  in  order  to  satisfy  the  circuit 
relation.  This  coronal  region  is  not  spatially  extended  because  the  pressure 
balance  assumed  in  its  construction  does  not  require  a  broad  spatial  region 
at  the  temperatures,  currents,  and  core  skin  depths  encountered.  The 
problem  of  experimental  examination  of  this  region  thus  becomes  more  deli¬ 
cate. 

For  the  model  as  presently  implemented,  it  is  possible  to  search  for 
surfaces  in  the  initial-condition  space  that  imply  the  bifurcation  of 
trajectories  among  the  SC/PC/8C  classes.  For  the  previous  version  of  the 
code,  this  was  not  useful  and  not  easily  done  because  of  the  extreme 
sensitivity  to  the  radiative  losses  discussed  above  in  regard  to  the  "exploded 
corona"  problem.  Such  partitioning  of  the  initial-condition  space  will 
provide  a  useful  guide  in  constructing  scaling  laws  and  in  estimating  the 
performance  of  alternate  generator  and  load  configurations,  because  the 
most  physically  appropriate  trajectories  (for  any  radiative  loss  approxima¬ 
tion)  will  be  those  near  the  PC/BC  bifurcation.  Moreover,  it  appears  that 
some  function  of  (f,  B,  r^^)^  will  provide  this  surface  (for  fixed  generator 
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parameters)  so  that  inference  back  to  initial  conditions  may  be  drawn 
from  later  performance,  or  the  domain  of  initial  conditions  limited  based 
upon  knowledge  of  particular  performance  criteria. 

In  addition  to  improved  radiative  loss  calculations,  there  remain  many 
areas  for  improvement  in  the  self-similar  model,  mo't  notably  the  inclusion 
of  magnetic  diffusion,  improved  dynamics  in  the  current-carrying  layer,  an 
equation  of  state  that  is  coupled  to  the  photoionization  process,  and  a 
continuous  radial  velocity  profile  for  the  core/corona  interface.  In  addi¬ 
tion,  the  performance  of  the  ID  code  suggests  that  some  direct 
dissipative  mechanism  for  radial  momentum  may  be  required  in  order  to  bring 
the  model  into  detailed  agreement  with  experiment.  If  such  a  mechanism 
is  found  then  the  PC/BC  bifurcation  surface  will  tend  to  widen  into  a  zone 
of  finite  measure  and  soften  the  present  sensitivity  to  the  initial  condi¬ 
tions.  These  problems  will  be  addressed  in  the  coming  months,  as  well  as 
extensions  of  the  core  self-similarity  model  to  allow  axial  or  axial/ 
azimuthal  variations. 
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Side  by  side  with  the  more  complete  analysis  embodied  in  the  SPLAT 
code  and  its  equations,  we  have  attempted  to  formulate  a  much  simpler  set 
of  model  equations,  simple  enough  to  give  various  analytic  scaling  laws 
for  the  collapse  time,  radiation  pulse  duration,  radiated  energy,  and 
temperature.  These  models,  necessarily  oversimplified  in  order  to  give 
analytic  results,  may  be  useful  for  making  some  further  approximations  in 
the  code,  and  for  predicting  order-of-magnitude  results  and  large-scale 
dependence  on  input  parameters  (mass  m,  initial  radius  r^,  atomic  weight  A) 
and  generator  parameters  (peak  voltage  V^,  voltage  rise  time  x^,  generator 
impedance  Z^). 

The  physics  of  the  system  depends  more  simply  on  the  current  I(t), 
peak  core  density  ni(t),  and  peak  implosion  velocity  ,  corona  'boundary' 
density  n,(t)  and  corona  scale  height  o(t)  _  all  of  which  are  related  to  the 

a 

input  and  generator  parameters  by  differential  equations  without  simple 
analytic  solutions.  It  is  our  hope  eventually  to  have  heuristic  fitting 
functions  relating  these  physical  parameters  to  the  input  and  generator 
parameters,  as  a  result  of  large  numbers  of  code  runs.  Difficulties  in  the 
code  development  have  prevented  our  having  enough  runs  yet  to  report  real 
achievement  of  this  objective,  although  we  do  not  see  fundamental  obstacles. 
Preliminary  scaling,  demonstrated  from  a  relatively  few  runs  is  shown  in 
Fig.  6. 

1 .  Collapse-Phase  Models. 

A  simple  model  represents  an  extreme  in  bounding  the  collapse- 
phase  behavior.  Plasma  pressure  and  corona  current  are  ignored,  and  the 
core-plasma  shell  is  accelerated  inward  by  only  the  ^xB^  force  due  to  current 
in  the  skin  layer,  whose  thickness  is  assumed  to  be  "frozen-in"  by  the  plasma 
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temperature  so  that  the  thickness  a(t)-r^(t)  =  f  a(t)  with  f  a  constant. 
The  equation  for  the  outer  radius  a  is  then 


-  a^  a  =  (TTmcf)“^  I^(t)  , 


if  plasma  compression  during  run-in  is  neglected  (until  assembly  on  the 
axis).  In  the  preliminary  code  runs,  I(t)  appears  to  be  nearly  linear  during 
the  collapse  phase,  i.e.,  the  effective  total  impedance  appears  to  be 
nearly  constant.  Using  the  simple  form  I  =  it  in  Eq.  (87)  gives  an  equation 
of  the  form  x^x  =  -  t^.  The  dimensionless  time  to  collapse  to  the  origin 
is  then  a  number  of  order  unity,  obtainable  by  numerical  solution,  and  the 
real  collapse  time  scales  as  a^I  ^(m  f)'‘.  The  more  complex  form  with  I  given 


Zgl  +  Lgl  +  (R  +  L)I  =  Vg(t) 

and  with  load  inductance  modeled  by  L  =  -  k  !ln(r/r^),  does  not  give 

readily  scalable  results. 


2,  Simplified  Assembly  Model. 

In  the  simplest  approximation,  we  may  imagine  a  more-or-less  uniform 
core  plasma  with  radius  r(t)  and  temperature  (T^  =  T^)  T(t),  with  the  origin 
of  t  taken  to  be  the  beginning  of  the  assembly  (i.e.,  we  assume  a  sharp 
assembly  beginning  at  t  =  o).  We  will  use  two  simultaneous  ordinary 
differential  equations  to  describe,  in  a  simplified  way,  the  evolution  of 
r(t)  and  T(t).  The  core  current  fraction  g  will  enter  in  general  as  a  para¬ 
meter  in  both  equations;  but,  to  the  extent  that  the  corona  has  robbed 
the  core  of  both  its  JxB.  compression  force  and  its  ohmic  heating,  the 
magnitude  and  effect  of  S  will  be  reduced. 
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The  force  equation  has  two  force  terms:  the  outward  2TTrp  =  2TrmT  with 


n  =  n.  +  Z(T)ng,  and  the  crrS^/Sir  magnetic  compression  force  with 
B  =  2SI{t)/cr(t) .  The  density  n-  is  given  by  N/TTr^(t)  with  fixed  N. 

The  energy  equation  has  an  assembly  heating  term  representing  the  con- 
version  of  kinetic  energy  H  mr^  into  thermal  energy  j(n^.  +  Zn^Tland  ioniza¬ 
tion  potential  energy  £^(2)n.T,  and  has  radiative  cooling  terms  corresponding 
to  line  radiation  with  a  line-average  opacity  factor,  recombination  radiation, 
and  bremsstrahlung,  which  will  be  allowed  to  escape  unattenuated.  The 
radiative  power  expressions  are  modeled  by: 

P^(W)  =  irr^  £-4.73  x  10‘^^n?Z  Fj_(T)  (88) 

-1  -U 

with  p^  =  .286  T  (£n  2t)  the  line  photon  escape  factor, 

'  “  m(g)/il(cm)]  , 


Fl(T)  approximating  the  results  of  Duston  &  Davis  (Ref.  19)  and 


~^keV  1  ^ 

(91) 

and  Z(t)arT|[/J  for  .03  <  .25  , 

(92) 

Pr(W)  =  Pg  x  Min  (3,  ZT'gJ/^) 

(93) 

Pg(W)  =  TTrH  X  4.74  X  10"^^  "i^’''’keV  ’ 

(94) 

Upon  specializing  to  Aluminum  (i;a!21)  and  expressing  all  times  in  'shakes', 
-8 

i.e.,  units  of  10  sec  (the  natural  units  for  the  problem),  we  have  the 
following  pair  of  simultaneous  ordinary  differential  equations: 


(89) 

(90) 


no 


* 


(95) 


r  =  Kpr-'' 

f  =-K^r  r  -  [k^  -  iin{r  T*'*)f  F^d)!''^  (96) 

-  K^r'^T^  -  Kjr-2T3/2  +  •  4  ^  d'' 

with  T  in  keV,  r  in  trin,  I  in  f^,  V  in  MV,  and  with: 

^5.5  (Pressure  expansion  factor)  (97) 

V  =  '-o  s'lj,  (Magnetic  compression  coef.)  (98) 

^  (Assembly  heating  coef,)  (99) 

■  1-97  X  10  ^  (Line  rad.  coef.)  (TOO) 

<3  =  10.9  d  ni/£  =  0.48u  (Bremsstrahlung  coef.)  (101) 

Kr  a  3  Kg  for  T  <  .3,  2  Kg  for  T  ^  .3  (Recomb.  rad.  (102) 

coef. ) 

e/10‘‘m(g)  (Ohmic  heating  coef.)  (103) 

’^T  "  +  in  u  (Optical  depth  factor)  (104) 

y  *  10'*m(g)/x,(cm)  (Mass  coef.)  (105) 

R  ~  3/2  for  T  <  .3,  7/6  for  T  >  .3  (Exponent  for  recomb.)  (106) 
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Once  the  temperatures  reach  0.24  keV,  the  and  f  are  replaced  by 

.62T  and  .627  respectively.  (This  assumes  that  all  the  atoms  are  in  the 

average  ionization  state  I  and  that  above  240  eV  the  A1  is  stripped). 

Here  the  coronal  effects  enter  via  the  unspecified  parameter  3=1  / 

core 

^total . 

From  the  coefficients  in  this  form,  one  sees  that  for  m-lO^'^g, 

£~3  cm,  '"i  ~  1  mm,  Tg~0.1  keV,  I~  1  MA  and  V~1  MV  one  will  have  "initially" 
at  the  beginning  of  assembly 


r  ~  [0.7  -  (3)6^]  tnm/sh^ 

-  k^r  r  ~  .07j  r-i  l(mm/sh)  •  (1  -  4S^) 

yv  ~  0.7S 

7^^^  ~  .005 

~  .015 

KLr'V|_(T)  T=  .023 

If  f  remained  constant,  a  minimum  radius 

'■n.in-'-l  -'i'-j  Co-?  - 

would  be  reached  in  a  time 


ta  ~  1^11:0.7  -  (3)6^]-'' 


acceleration 

PdV  heating,  if  ^r‘)<0 
Ohmic  heating 
Bremsstrahlung  cooling 

Recomb.  rad.  cooling 

Line  rad.  cooling 


(107) 


(108) 
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If  radiation  cools  T  faster  than  r"^  can  increase,  then  the  is  smaller 
(cooling  then  even  faster)  and  t.  is  shorter;  so  the  condition  for  refrigera- 
tive  collapse  with  small  S  is  approximately  that 

K„r“^  decreases  on  the  timescale  IrJ/K  r"^ 

P  ‘  r  p 


i.e. , 


ii  +  Liril 

3  T  r  1 


>  I  K  r-1 
3  p 


(109) 


assuming  !4f/3T|  >  (-r)/r.  This  is  clearly  easier  when  [r^|  is  large. 
When  such  a  collapse  occurs,  all  the  kinetic  energy  of  implosion  is  con 
verted  into  radiation,  as  well  as  the  initial  thermal  and  ionization 
energy.  The  fate  of  the  magnetic  energy  depends  on  the  behavior  of  3l 
during  such  a  collapse. 


But  the  amount  of  radiated  energy  coming  out  at  energies  above,  say 
1  keV,  depends  mostly  on  the  history  of  radiation  before  the  plasma  cools 
to  T~300  eV,  and  the  exponentially  small  number  of  ^  1  keV  transitions 
excited  in  a  plasma  of  T  <  300  eV.  Thus,  the  experiments  optimizing 
radiation  above  1  keV  will  be  those  that  get  the  plasma  hotter  than  about 
300  eV  before  or  during  collapse,  and  have  low  enough  densities  at 
assembly  to  prevent  radiation  cooling  below  300  eV  on  a  fast  timescale. 
Regarding  this  threshold  temperature  T^  as  a  parameter,  one  wishes  to 
maximize 

fit  [K|_  r'*  r-'  , 

CT  >  T,: 

in  other  words,  if  S  is  small,  maximize 

CT  >  Tg] 
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Here  -f/r  is  approximately  the  inverse  of  the  compression  timescale,  and 
-f/T  the  inverse  of  the  net  cooling  timescale  (dominated  by  radiative 
cooling). 

Note  that,  because  of  the  large  opacity  factor  in  the  line  radiation 
and  the  V  dependence  of  bremsstrahlung  (as  opposed  to  dependence  in  the 
line  radiation),  the  line  and  continuum  cooling  rates  may  be  comparable. 

Note  also  that  the  simultaneous  equations  for  r  and  T  in  even  this  over¬ 
simplified  model  are  complex  enough,  and  F|^(T)  is  complex  enough,  that  no 
simple  scalings  beyond  the  ones  just  discussed  can  be  obtained.  As  a 
simple  rule,  we  have  a  radiation  pulse  length 

ta-lril/r-lr^lr^  Kp’  (m) 

and  a  total  radiated  energy  of  order  h  nirj,  but  the  radiative  yield  above 

a  threshold  T„  cannot  be  derived  as  a  simple  function  of  m,  J,,  V,  I,  and 
y 

the  generator  risetime. 

The  differential  equations  of  this  model  were  solved  for  certain  cases 
with  only  bremsstrahlung  cooling  and  with  8  =  0.  "Initial"  velocities  just 
prior  to  collapse  were  set  at  values  of  order  10^  cm/sec  and  initial 
temperatures  at  values  of  order  0.1  keV.  Transparent  bremsstrahlung  alone 
was  then  enough  to  induce  simple  collapse  trajectories,  in  which  radiation 
cooling  dominated  the  thermal  energy  input  from  the  collapse.  The  decelera¬ 
tion  due  to  plasma  pressure  even  without  any  opposing  ^xB  force  was 
insufficient  (by  an  order  of  magnitude)  to  prevent  the  density  buildup 
responsible  for  the  cooling.  In  order  to  do  so,  the  peak  density-squared, 
proportional  to  -  (~l)rp"^»  had  to  be  lowered  by  two  orders  of 
magnitude  (as  achievable  when  r^~0.3  mm/shake). 
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3.  Radiation  Scaling  Model. 

Mosher  (Ref.  5)  has  derived  radiation  scaling  with  wire  mass  m  and 
assembled  radius  ri ,  based  on  the  following  assumptions: 


(1)  Energy  input  is  the  implosion  kinetic  energy  K  =  ’5  mv|, 
delivered  to  the  plasma  on  the  hydrodynamic  assembly 
timescale  t^^  = 

(2)  Energy  loss  is  by  radiation,  to  which  the  plasma  is 
transparent  above  temperature  T^.  The  radiated  power 
is  modeled  by 


^rad  "  '’bb  "  '^BB’'i''’'‘  ^  ’’’c 


(112) 


2  .-2  .-0.7 


Prad  =  ^  T 


for  T  >  T, 


Tj.  is  the  value  at  which  the  two  expressions  are  equal. 


(3)  The  heating  power,  ^  exceeds  the  maximum 


so  that  T  rises  to  exceed  T^. 


(4)  Energy  not  radiated  is  stored  in  internal  energy  and  in 
temperature:  e  =  Kg.mT,  with  c  =  P^^  -  P^gj-  The  model, 
however,  appears  to  give  a  yield  scaling  with  m,  r.j ,  and 
v.|  which  disagrees  with  experiments.  The  model  may  be  mis 
leading  in  that 

(a)  it  neglects  ohmic  heating,  which  continues  after 
assembly  and  could  account  for  most  of  the  energy 
input  after  assembly. 
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(b)  there  is  evidence  that  the  energy  loss  precedes  trans¬ 
parency  of  the  plasma;  in  fact,  the  plasma  is  probably 
never  transparent,  but  nor  is  it  a  black  body  for  the 
dominant  lines. 

(c)  Continuum  radiation  is  neglected  in  the  T  >  domain, 

(d)  no  account  is  taken  of  coronal  effects,  and 

(e)  the  single  differential  equation  for  T  which  results, 

K_m  t  =  I  <  t^)  - 

assumes  a  constant-density  radiator,  i.e.,  the  radiation 
cooling  rate  cannot  induce  contraction. 

We  have  modified  this  model  to  include  ohmic  heating,  continuum 
radiation,  and  a  'typical'  opacity  for  the  line  radiation.  Coronal 
effects  can  be  studied  in  part,  through  the  parameter  3  in  Eqs.  98  and  103. 
Eqs.  95  -  106  remove  the  disadvantages  just  mentioned,  but  do  not  give 
particularly  tractable  results  useful  for  scaling  laws.  By  taking  the 
assembly  heating  to  be  'imv^/t^  with  t^  =  ir-jj/Er,  and  eliminating  r,  a 
generalization  of  Eq.  (113)  can  be  derived  from  (95)  -  (106): 

=  KpK^r'^  -  Kj_K^r'''  -  3Kgr'^  +  Kj^VI(t)  (114) 

with  Fj^(T)  approximating  the  results  of  Duston  and  Davis  (Ref.  19).  Now 

the  assembled  radius  r  is  simply  a  parameter,  as  in  the  Mosher  analysis, 

and  not  a  time-dependent  variable  with  an  equation  governing  its  evolution. 

At  the  temperatures  expected,  F|^(T)  is  non-monotonic  but  F^^(T)T'^ 

-2  3/2 

+  3  Kgr  T  may  well  be.  Approximating  the  assembly  and  radiation 
T-dependence  by  a  single  power  law,  one  gets 
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(115) 


TssK^T  +  Kj^VI(T)r’^^ 

with  KpK^r"^  T  -  F|^(T)T^''®  -  3  Kgr'^T^'^®  modeled  by  K^T.  The 

time-dependence  of  I  is  ignored,  and  I  is  taken  as  a  parameter  of  the 
problem  (hiding  its  T  dependence).  The  solution  of  Eq.  (115)  for  k.|>0  is 

-  K^/'Kj^  ;  (116) 


r4/3  _ 


K2/K., 


+  T 


4/3 


exp 


with  <2  =  KoVI.  For  K.|  <  0, 

=  K2/!,<t!  -  [k2/!K^|  -  T^^^jexp  (-  f  .'K^'t) 


(117) 


The  reference  temperature  T.j  may  be  taken  small  and  neglected.  Doing  so 
gives 


T  =  (1^2/ 1  I 


(118) 


Ordinarily  is  negative,  and  Eq.  (118)  describes  the  Ohmic  heating  to 

final  temperature  (K2/!K.jj)  in  the  presence  of  radiative  losses,  until 
the  assembly  heating  ceases.  Then  with  T  as  a  new  "initial"  temperature 
T-],  Eq.  (117)  with  the  new  K-j  describes  the  subsequent  temperature  decay. 
The  central  parameter,  K2/lK^),  scales  as  SlVm“^|K^l’^  with 

|K^|  =  K^_K  (T^)r'^  Fj_(T^)T^''^  +  ^K^r'^Tj^®  -  KpK^r"^  (119) 


and  the  remaining  coefficients  given  by  Eqs.  (97)  -  (106).  From  this 
prescription,  the  duration  for  which  T  >  Tg  (some  threshold  temperature] 
can  be  calculated,  and  the  radiation  yield  during  that  time  interval  can 
be  calculated  from  T(t).  The  assembled  plasma  radius  remains  as  a  free 
parameter,  which  is  the  price  paid  for  ignoring  its  simultaneous  evolution. 
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FIGURE  CAPTIONS 


Figure  1  -  Evolution  of  the  Iniplicit-T2  Model.  The  inner,  outer  core 

annulus  boundaries  Cb,a]  are  in  (cm),  all  other  quantities  are 
in  arbitrary  units.  The  time  is  in  units  of  the  assembly  time, 
t  :  b  =  0. 

a 

Figure  2a-  Annulus  Evolution  and  Radiated  Power  in  the  Explicit-T2  Model. 
The  boundaries  [b,al  are  in  (cm),  is  normalized  to  its 
maximum. 

Figure  2b-  Thermodynamic  Variables  and  Circuit  Response.  All  quantities 
are  in  arbitrary  units. 

Figure  3a-  Two  "initial"  ion  density  profiles  for  different  SPLAT-ID 
3b 

code  runs.  Initial  temperatures  were  taken  to  be  15  eV. 

The  assumed  initial  current-penetration  radius,  r^ ,  is  shown 
for  both. 

Figure  -la-  Simple  collapse  (SC)  evolution  of  r^ ,  t^^,  b,  I,  and  Ppgjj- 
Insert  shows  the  relation  of  P^^^  to  ion  temperature  and 
spatially  averaged  ion  density. 

Figure  4b-  Pause-and-col lapse  (PC)  evolution  of  r^. ,  r^,  b,  I,  and 
Insert  shows  the  relation  of  P^^^j  to  ion  temperature  and 
spatially  averaged  ion  density. 

Figure  5  -  Generator  voltage  vs  time  for  a  =  1.2  (1.2  x  faster  than  normal 
pulse) . 
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gure  6  -  Frequency-integrated  radiative  yield  vs  generator  rise-rate 
multiplier  a,  for  three  initial  configurations:  (a)  that  of 
Fig.  3a,  (b)  that  of  Fig.  3b,  and  (c)  a  smaller  radius  array 
with  peak  at  0.3  cm  and  width  0.1  cm,  but  the  same  mass  as 
3a  and  3b. 


122 


DISTRIBUTIOfi  LIST 


3EPARTMENT  OF  OEFENSE  DEPARTMENT  OF  THE  NAVY  (Continued) 


Assistant  to  tne  Secretary  of  Defense 

Naval  weapons  Ctr 

Atomic  Energy 

ATTN 

Code  343,  FKA6A2,  Tech  Svcs 

ATTN 

Executive  Assistant 

ATTN 

Military  Applications 

DEPARTMENT 

OF  THE  AIR  FORCE 

Defense  Intelligence  Agency 

Air  Force  Heaoons  Lao 

ATTN 

DT-IB,  R.  Rubenstein 

ATTN 

SUL 

ATTN 

CA 

Defense  Nuclear  Agency 

ATTN 

STYC 

ATTN 

RAAE 

ATTN 

NT 

ATTN 

RAEV 

ATTN 

\TYP 

ATTN 

STNA 

ATTN 

RAEE 

Ball istic  Missile  Ofc 

4  cy  ATTN 

TITL 

ATTN 

ENSN 

ATTN 

SYOT 

Defense  Teen  Info  Ctr 

ID  cy  ATTN 

DO 

Deputy  Chief  of  Staff 

Rsch,  Dev, 

4  Acq 

'ieli  Command  Defense  Nuclear  Agency 

ATTN 

AFROQl 

3et  1 

Dawrence  Livermore  LaO 

Space  Oiv 

ATTN 

FC-1 

AUN 

XR,  Plans 

ATTN 

YEZ 

Field  Command 

ATTN 

YGJ 

Defense  Nuclear  Agency 

ATTN 

YKF 

AT'N 

FCPR 

ATTN 

YKS,  P.  Stabler 

ATTN 

FCTT,  W,  Sunina 

ATTN 

YKM 

ATTN 

FCTT 

ATTN 

YNV 

ATTN 

FCTXE 

ATTN 

FCTT,  G.  Ganong 

Strategic  Air  Command 

ATTN 

INT.  E.  Jacobsen 

Jnder  Secretary  of  Defense  for  Rsen  4  Engrg 

ATTN 

strategic  4  Space  Sys  (OS) 

DEPARTMENT  OF  ENERGY 

DEPARTMENT 

OF  the  army 

Department  of  Energy 

Office  of  Military  Application 

Harry  Diamond  Laos 

ATTN 

Ofc  of  Inert  Fusion,  S.  Kanalas 

ATTN 

OELHO-NW-P,  20240 

ATTN 

Ofc  of  Inert  Fusion,  T.  GoQlove 

ATTN 

OELHO-NW-RA,  22100 

ATTN 

Ofc  of  Inert  Fusion,  C.  riiHand 

ATTN 

OELHO-.NW-RI,  Kervis,  22900 

ATTN 

OELHO-TA-L.  31100 

OTHER  GOVERNMENT  AGENCY 

JS  Army  Nuclear  4  Chemical  Agency 

Central  Intelligence  Agency 

ATTN 

Li brary 

ATTN 

0SHR7NED 

JS  Army  Test  and  Evaluation  Comd 

DEPARTMENT  OF  ENERGY  CONTRACTORS 

ATTN 

ORSTE-EL 

University 

Of  California 

JSA  Missile  Command 

Lawrence  Livermore  National  Lab 

ATTN 

Ooc  Section 

ATTN 

Tech  Info  Dept,  Library 

AUN 

L-S45,  J.  Nuckolls 

DEPARTMENT 

OF  the  navy 

ATTN 

L>47,  L.  Uouters 

ATTN 

L-153 

Naval  Rsch 

Lab 

ATTN 

L-153,  0.  Meeker 

ATTN 

Code  2000,  J.  Brown 

ATTN 

M.  Pickles.  L-401 

ATTN 

Code  4770,  I.  Vltokovitsky 

ATTN 

Code  4720,  J.  Oavis 

Sandia  National  Lab 

ATTN 

Code  4780,  S.  Ossakow 

ATTN 

G.  Yanas 

ATTN 

Code  4773,  G.  Cooperstein 

ATTN 

Tech  Lib,  3141 

ATTN 

J.  Powell 

Naval  Surface  Weapons  Ctr 

ATTN 

N.  Clauser,  Org  5241 

ATTN 

Code  R40 

ATTN 

G.  Kuswa,  Org  5240 

AHN 

Code  F31 

AHH 

0.  Allen,  Org  9335 

ATTN 

Code  F34 

123 


DEPARTMENT  OF  ENERGY  CONTRACTORS  IContinuea; 

Los  Alamos  National  Lao 

ATTN:  MS222,  J.  Srownell 

DEPARTMENT  OF  DEFENSE  CONTRACTORS 

Advanced  Researcn  i  Applications  Corp 


ATTN 

R.  Armistead 

Aerospace  Corp 

ATTN 

V.  jQsepnson 

ATTN 

S .  Bower 

ATTN 

*ecn  Info  Svcs 

80H  Corp 

ATTN 

Corporate  LiPrary 

SOM  Corp 

ATTN 

L.  Hoeft 

Boeing  Co 

ATTN 

Aerospace  LiOrary 

Oi kewood 

ATTN 

Tech  Lib  for  L.  Davis 

EG&G  Wash  Analytical  Svcs  Ctr,  Inc 

ATTN 

Li  brary 

General  Electric  Co 

ATTN 

J.  Peden 

ATTN 

H.  O'Donnell 

IRT  Corp 

ATTN 

R.  Mertz 

JAVCOR 

ATTN 

E.  wenaas 

JAYCOR 

ATTN 

E.  Alcaraz 

ATTN 

R.  Sul  1 ivan 

A  cy  ATTN 

J.  Guillory 

A  :y  ATTN 

R.  Terry 

<aman  Sciences  Corp 
ATTN:  S.  Face 


<aman  Tempo 

ATTN:  DASIAC 


DEPARTMENT  OF  DEFENSE  CONTRACTORS  'Continued, 

Lockneed  Missiles  i  Space  Co,  Inc 
ATTN:  L.  Chase 

Lockneed  Missiles  1  Space  Co,  Inc 

ATTN:  5.  Taimuty,  Dept  31-7A/1SA 


Maxwell  Laps,  Inc 


ATTN 

A  .  <0  1  b 

ATTN 

3.  'anmoto 

ATTN 

A.  Miller 

ATTN 

■) .  Co  ‘  e 

'ic^onne!  1 

'ouglas  Cor? 

ATTN 

S.  Scnneider 

Mission  Asch  Corp 

ATTN 

W.  Hart 

ATTN 

C.  Longmire 

Mission  Rscn  Corp 

ATTN 

3.  jOOfrey 

Mission  Rscn  Carp,  San  Diego 

ATTN 

3.  Passenheim 

Pacific-Sierra  Rscn  Cor? 

ATTN 

H.  Brode,  Chairman 

ATTN 

L.  Schlessinger 

^nvsics  International  Co 

ATTN 

C.  jUman 

ATTN 

C.  Stal 1 ings 

ATTN 

S.  "razier 

R  4  D  Associates 

ATTN 

HJas 

ATTN 

A.  Latter 

R  4  D  Associates 

ATTN 

P.  Turcni 

S-CUBED 

ATTN 

A.  Wilson 

Sciences  Applications,  Inc 

ATTN 

K.  Sites 

TRW  Electronics  4  Defense  Sector 
ATTN:  Tech  Info  Ctr 
ATTN:  D.  Clement 


124 


